NPS~63Cj 78031. 


NAVAL  POSTGRADUATE  SCHOOL 


// 


Monterey,  California 


FINAL  TECHNICAL  REPORT  OF  NSF  GRANT  DES  75-10719 
"TROPICAL  WAVE  DYNAMICS" 

by 

Chlh-Pei  Chang  and  Roger  T.  Williams 
March  1978 

Final  Report  for  Period     October  1975  -  March  1978 


FEDDOCS 

D  208.14/2:NPS-63CJ78031 


Approved  for  public  release;  distribution  unlimited, 
spared  for:   National  Science  Foundation 


Washington,  D.C.   20550 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 

Rear  Admiral  Tyler  F.  Dedman  J.  R.  Borstin^ 

Superintendent  Provost 


The  work  reported  herein  was  supported  by  the  National  Science  Founda- 
tion, Atmospheric  Research  Section,  under  Grant  DES  75-10719. 

Reproduction  of  all  or  part  of  this  report  is  authorized. 

This  report  was  prepared  by: 


UNCLASSIFIED 


SECuBlTY    Cl  ossification   OF    THIS   PAOE   (Whon   Date  Entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


I.     REPORT   NUMBER 

NFS  63Cj 78031 


2.   GOVT   ACCESSION   NO 


3.     RECIPIENT'S  CATALOG   NUMBER 


4.     Title  rand  Subf/f/e) 

Final  Technical  Report  of  NSF  Grant  DES  75-10719 
"Tropical  Wave  Dynamics" 


5.     TYPE  OF   REPORT   ft   PERIOD  COVERED 

Final  Report  for 
Oct  1975  -  Mar  1978 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7   AUTHORCs; 

Chih-Pei  Chang  and  Roger  T.  Williams 


8.  CONTRACT  OR  GRANT  NUMaERfaJ 

NSF  DES  75-10719 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  California   93940 


10.     PROGRAM   ELEMENT.  PROJECT,   TASK 
AREA  &   WORK  UNIT  NUMBERS 


11.     CONTROLLING  OFFICE  NAME   AND   ADDRESS 

National  Science  Foundation 
Washington,  D.C.   20550 


12.     REPORT  DATE 

March   1978 


13.     NUMBER  OF  PAGES 

65 


14.     MONITORING  AGENCY  NAME  4   ADDRESSC// d///er»nf  Irom  Controlling  Ofllco) 


15.     SECURITY  CLASS,  (ol  Ihia  roport) 

Unclassified 


I5«.     DECLASSIFI  CATION/ DOWN  GRADING 
SCHEDULE 


16.     DISTRIBUTION   STATEMENT  (ol  thia  Report) 

Approved  for  public  release;  distribution  unlimited. 


17.     DISTRIBUTION  STATEMENT  (ol  the  abatfct  ont»r»d  In  Block  30,   II  dllltrent  Irom  Raporl) 


18.     SUPPLEMENTARY  NOTES 


19.     KEY  WORDS  (Contlnuo  on  revarae  aide  II  necaaattry  and  Idontify  by  block  numbmr) 

Tropical  Waves 

Tropical  Meteorology 

Tropical  Dynamics 

Monsoon 

Barotropic    Instability 


20.     ABSTRACT  (Contlnua  on  ravaraa  aldm  it  nacaaaary  and  Identify  by  block  number) 

This  research  may  be  divided  into  the  following  three  parts:   a)  waves 
forced  by  condensation  heating  in  the  troposphere,  b)  the  local  barotropic 
instability  produced  by  a  spatially  varying  basic  flow,  c)  the  numerical 
modeling  of  the  wave-basic  flow  interacting  system.   In  the  first  part,  the 
problem  of  the  vertical  structure  of  tropical  waves  generated  by  cumulus 
heating  is  treated  for  the  following  two  cases:   1)  heating  is  controlled  by 
large-scale  low-level  convergence  and  2)  heating  receives  no  feedback  from 


00     )    1AM   73     1473  EDITION  OF    1  NOV  65  IS  OBSOLETE 


S/N    0  102-0J4- 660  1 


UNCLASSIFIED 


SECUBITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  Bnlarad) 


UNCLASSIFIED 


ii^lTY    CLASSIFICATION  OF   THIS  P  XQEfWhan  Data  Entered) 


large-scale  motion.   In  the  first  case  it  is  shown  that  waves  which  have  a 
short  vertical  wavelength  relative  to  the  vertical  scale  of  heating  are 
stable,  and  the  structure  of  the  most  unstable  solution  is  computed.   In  th 
second  case  it  is  shown  that  the  forcing  is  most  efficient  for  the  longest 
zonal  wavelength  and  that  the  most  favored  vertical  wavelength  is  about  twic 
the  vertical  scale  of  the  heating.   The  effect  of  cumulus  damping  on  plane- 
tary-scale Kelvin-like  oscillation  gives  solution  which  resembles  certain 
observed  tropical  circulation. 

The  second  portion  of  the  research  treats  waves  in  a  barotropically  un- 
stable jet,  which  varies  in  x  and  y.   Numerical  solutions  are  obtained  in 
a  channel  with  time  periodic  forcing  on  the  east  side  and  a  radiation  condi 
tion  on  the  west  side.   The  solutions,  which  are  periodic  in  time  with  the 
forcing  frequency,  show  the  spatial  growth  and  decay  of  the  waves  as  they 
move  through  the  channel. 

In  the  third  portion  of  the  research  a  3-level  primitive  equation  model 
is  used  to  simulate  the  interaction  between  the  planetary-scale  monsoon  cir 
culation  and  the  synoptic  waves.  The  heating  is  specified  from  the  200-mb 
time-mean  divergence  field  analyzed  by  Krishnamurti  and  Rogers.  The  inter- 
actions are  examined  with  a  constant  heating  and  with  a  heating  with  a  ten 
day  period.  The  model  properly  simulates  certain  aspects  of  the  summer  mon 
soon  circulation  but  also  reveals  certain  problems  in  the  vorticity  damping 
mechanism. 


UNCLASSIFIED 


SECURITY  CLASSIFICATIOU  OF  THIS  PAGE(TW««n  Datm  Bnlarad) 


NATIONAL  SCIENCE  FOUNDATION 
Washington,  D.C.    20550 


SUMMARY  OF  COMPLETED  PROJECT 


Form  Approved 
OMBNo.  99R0013 


Please  read  instructions  on  reverse  carefully  before  completing  this  form. 


1.  INSTITUTION  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  Calif  939A0 


2.  NSF  PROGRAM 

Meteorology 


3.    PRINCIPAL  INVESTIGATOR (S) 

Chih-Pei  Chang 
Roger  T.  Williams 


4.    AWARD  NUMBER 

DES   75-10719 


5.    DURATION 

.MOS,  3Q 


6.    AWARD  PERIOD 

fromlO/lpSto    3/31/78 


7.    AWARDEE  ACCOUNT  NUMBER 

56522 


&    PROJECT  TITLE 


Tropical  Wave  Dynamics 


9,    SUMMARY  (ATTACH  LIST  OF  PUBLICATIONS  TO  FORM) 

This  research  may  be  divided  into  the  following  three  parts:   a) waves 
forced  by  condensation  heating  in  the  troposphere,  b)  the  local  barotropic 

istability  produced  by  a  spatially  varying  basic  flow,  c)  the  numerical 
modeling  of  the  wave-basic  flow  interacting  system.   In  the  first  part,  the 
problem  of  the  vertical  structure  of  tropical  waves  generated  by  cumulus 
heating  is  treated  for  the  following  two  cases:   1) heating  is  controlled  by 
large-scale  low-level  convergence  and  2)  heating  receives  no  feedback  from 
large-scale  motion.   In  the  first  case  it  is  shown  that  waves  which  have  a 
short  vertical  wavelength  relative  to  the  vertical  scale  of  heating  are 
stable,  and  the  structure  of the  most  unstable  solution  is  computed.   In  the 
second  case  it  is  shown  that  the  forcing  is  most  efficient  for  the  longest 
zonal  wavelength  and  that  the  most  favored  vertical  wavelength  is  about  twice 
the  vertical  scale  of  the  heating.   The  effect  of  cumulus  damping  on  plane- 
tary-scale Kelvin-like  oscillation  gives  solution  which  resemble  certain  ob- 
served tropical  circulations. 

The  second  portion  of  the  research  treats  waves  in  a  barotropically  un- 
stable jet,  which  varies  in  x  and  y.  Numerical  solutions  are  obtained  in 
a  channel  with  time  periodic  forcing  on  the  east  side  and  a  radiation  condi- 
tion on  the  west  side.   The  solutions,  which  are  periodic  in  time  with  the 
forcing  frequency,  show  the  spatial  growth  and  decay  of  the  waves  as  they 
move  through  the  channel. 

In  the  third  portion  of  the  research,  a  3-level  primitive  equation  model 
is  used  to  simulate  the  interaction  between  the  planetary-scale  monsoon  cir- 
culation and  the  synoptic  waves.   The  heating  is  specified  from  the  200-mb 
time-mean  divergence  field  analyzed  by  Krishnamurti  and  Rogers.   The  inter- 
actions are  examined  with  a  constant  heating  and  with  a  heating  with  a  ten 
day  period.   The  model  properly  simulates  certain  aspects  of  the  summer  mon- 
soon circulation  but  also  reveals  certain  problems  in  the  vorticity  damping 
mechanism. 


9.  SIGNATURE  OF  PRINCIPAL  INVESTIGATORS 
PROJECT  DIRECTOR 


TYPED  OR  PRINTED  NAME 

Chih-Pei  Chang 


DATE 


NSF  Form  98A  (10-77)    SUPERSEDES  ALL  PREVIOUS  EDITIONS 


TABLE  OF  CONTENTS 

1.  Description  of  Research  and  Results  ------------   5 

2.  List  of  Publications   -------------------12 

3.  Theses 13 

4.  Scientific  collaborators   -----------------  13 

5.  Comments   -------------------------13 

6.  Appendix  (reprints)  --------------------14 

Initial  Distribution  List  -------------------65 


1.   DESCRIPTION  OF  RESEARCH  AND  RESULTS 

The  research  goal  of  this  project  may  be  summarized  into  the 
following  three  parts 

a)   Waves  forced  by  condensation  heating  in  the  troposphere 
We  have  studies  the  problem  of  the  vertical  structure  of  trop- 
ical waves  generated  by  cumulus  heating  for  two  cases:   l)heating 
is  controlled  by  large-scale  low-level  convergence  (CISK  model) 
and   2)  heating  receives  no  feedback  from  large-scale  motion 
[forced  model) .   In  the  first  case  the  internal  wave-CISK  system 
is  treated  and  upper  and  lower  bounds  for  the  vertical  wavelength 
of  the  unstable  waves  under  normal  heating  conditions  are  esta- 
blished through  analysis  of  the  frequency  equation.   The  lower 
bound  excludes  the  possibility  of  excitation  or  maintenance  of 
short  vertical  wavelengths  relative  to  the  vertical  scale  of  heat 
ing.   The  vertical  structure  of  the  most  unstable  waves  is  also 
computed.   These  results  were  reported  in  Chang  C1976a) .   In  the 
second  case  the  problem  of  scale-selection  of  Kelvin  waves  in  the 
stratosphere  is  analyzed  using  an  inviscid  linear  model.   The  re- 
sults modify  Holton's  (1973)  theory  in  that   1)  the  forcing  is 
most  efficient  for  the  longest  zonal  wavelength  even  if  the  heat 
sources  are  distributed  randomly,  and   2)  the  most  favored  verti- 
cal wavelength  of  the  excited  waves  is  about  twice  the  vertical 
scale  of  heating.   These  results  were  reported  in  Chang  C1976b). 
We  have  also  treated  the  linear  model  with  friction  added  which 
is  called  for  by  several  budget  studies  of  cumulus-related  tropi- 
cal weather  systems,  including  the  monsoon  CHolton  and  Colton, 
1972),  the  cloud  clusters  (Williams  and  Gray,  1973),  and  the 


easterly  waves  (Reed  and  Recker,  1971;  Wallace,  1971,  Reed  and 
Johnson,  1974;  Chu ,  1976).   This  friction  is  believed  to  be  due 
to  vertical  cumulus  transports.   Our  analysis  shows  that  this 
damping  has  a  strong  influence  at  low  frequencies  on  the  forced 
equatorial  waves  and  results  in  two  types  of  dispersive  relation- 
shj-ps .   The  first  type  is  characteristic  of  the  regular  internal 
gravity  waves  which  have  fast  phase  speeds  and  weak  vertical  atten-  I 
uation.   The  second  type  is  dominated  by  the  viscous  damping  time 
scale  and  has  slow  phase  speed  and  strong  vertical  trapping.   This 
theory  predicts  that  the  stratospheric  oscillations  may  be  identi- 
fied with  the  first  type  and  the  tropospheric  oscillations  with 
the  second  type.   In  the  case  of  Kelvin  waves  the  results  can  be 
used  to  explain  consistently  both  the  observed  stratospheric  Kelvin 
waves  and  the  planetary-scale  Kelvin-like  oscillations  in  the 
troposphere  including  the  40-50  day  oscillation  and  the  winter  mon- 
soon and  Walker  circulations.   These  results  are  described  in 
Chang  (19  77a). 

b)   The  local  barotropic  instability  produced  by  the  spatially 
varying  basic  flow. 

Here  we  have  studied  the  behavior  of  waves  superimposed  upon 
a  barotropically  unstable  mean  wind  field,  which  varies  in   x   and 
y.   This  mean  wind  field  roughly  simulates  the  200  mb  easterly  jet 
south  of  the  Tibetan  high  during  the  northern  hemisphere  summer. 
The  barotropic  vorticity  equation  is  linearized  and  a  linear  damp- 
ing term  is  added.   The  following  mean   u   component  is  used 

uCx,y)  =  -U(x)  sech^(y/d(x))  -  U^  , 


where   dCx)   is  a  measure  of  the  width  of  the  jet.   The  quantity 
U(x)  and  u   are  chosen  in  such  a  way  that  the  mean  flow  is  non- 
divergent.   A  rectangular  domain  is  employed  with  a  time  periodic 
boundary  contition  on  the  east  and  a  radiation  condition  on  the 
west,  which  simulates  the  propagation  of  small  amplitude  waves 
through  the  easterly  jet  region.   The  vorticity  equation  is  solved 
with  the  use  of  finite  differences,  and  when  the  boundary  condi- 
tions are  properly  adjusted  the  waves  move  smoothly  across  the  re- 
gion and  out  the  western  boundary.   After  a  certain  time  interval 
of  numerical  integration,  the  solution  becomes  periodic  everywhere, 
with  the  frequency  that  is  specified  on  the  eastern  boundary.   As 
the  waves  move  through  the  domain  they  grow  or  decay  spatially  in 
reaction  to  the  local  basic  flow,  even  though  at  each  point  the 
variation  is  purely  harmonic.   We  have  developed  a  simple  mechan- 
istic analytical  model  which  describes  the  main  features  of  the 
spatial  variation.   The  simple  model  uses  growth  rates  and  phase 
velocities  which  are  computed  from  a  parallel  flow  model  with  the 
wavelengths  taken  from  the  full  numerical  model.   These  wave- 
lengths increase  to  a  miximum  at  the  point  where  the  jet  is  most 
unstable,  and  then  decrease;  the  phase  velocity  has  a  similar  be- 
havior.  The  maximum  amrjlitude  in  the  numerical  solution  occurs 
downstream  from  the  most  unstable  part  of  the  jet,  at  the  point 
where  the  growth  rate  in  the  parallel  flow  model  becomes  zero. 
The  above  results  are  grossly  similar  to  that  expected  from  the 
parallel  flow  theory  of  barotropic  instability.   However,  in  the 
unstable  region  the  resultant  structure  of  the  waves  causes  a  spa- 
tial growtfi  rate  greater  than  predicted  by  the  local  growth  rates 


computed  with  a  parallel  flow  model.   In  the  stable  region,  the 
structure  leads  to  a  strong  dynamic  damping  which  is  apparently- 
due  to  the  presence  of  the  continuous  spectrum  modes  (Case,  1960; 
Pedlosky,  1964).   When  a  uniform  advective  velocity  is  added  to  a 
variable  mean  flow,  the  difference  between  the  behavior  of  the 
computed  waves  and  that  implied  by  the  parallel  flow  theory  is 
somewhat  reduced.   However,  in  this  case  a  stronger  zonal  asymmetry 
in  the  spatial  growth  rate  with  respect  to  the  jet  maximum  occurs 
as  a  result  of  slower  adjustment  of  the  wave  structure  to  the  local 
stability  conditions.   Chang  and  Williams  (1976)  presented  the 
analytic  model  and  discussed  some  other  analytical  aspects  of  the 
problem.   The  numerical  model  results  and  comparisons  with  the 
analytic  model  will  be  reported  by  Tupaz,  Williams  and  Chang  (.1978). 

c)   The  numerical  modeling  of  the  wave-basic  flow  interacting 
system. 

In  this  part  we  adopted  a  simple  numerical  model  formulated 
by  Monaco  and  Williams  (19  75)  to  simulate  the  planetary-scale  mon- 
soon circulation  and  the  interaction  with  synoptic  waves.   In 
order  to  concentrate  on  the  vorticity  dynamics,  the  heating  is 
specified  with  a  horizontal  distribution  given  by  the  time-mean 
200-mb  divergence  of  the  1967  summer  analyzed  by  Krishnamurti  and 
Rogers  (1970),   The  domain  of  the  model  extends  from  18°S   to 
46°N   and  there  are  three   a-levels  in  the  vertical.   Topography 
is  not  included  although  its  thermal  effect  is  implied  in  the 
heating  function.   Two  preliminary  experiments  have  been  carried 
out  to  study  the  vorticity  budget  of  the  planetary  scales  and  the 
non-linear  barotropic  energy  conversions.   In  one  experiment  the 


heating  is  fixed  in  time  while  in  the  other  it  fluctuates  with  a 
period  of  10  days.   The  time-mean  fields  with  moderate  damping  in 
both  experiments  show  a  significant  westward  phase  shift  of  the 
Tibetan  high,  which  suggests  that  the  correlation  of  the  fluctua- 
tions of  the  divergence  and  vorticity  as  suggested  by  Holton  and 
Colton  (1972)  is  probably  not  sufficient  for  the  required  planetary' 
scale  vorticity  sink.   However,  the  integration  period  is  only 
three  times  the  fluctuation  period  in  the  second  experiment  and 
the  amplitude  of  fluctuation  is  about  half  that  observed.   So  the 
full  effect  of  the  fluctuation  may  not  have  been  properly  repre- 
sented in  the  experiment.   A  longer  integration  with  larger  fluc- 
tuation is  probably  necessary  to  assess  fully  the  effect  of  the 
fluctuating  heating.   On  the  other  hand,  the  positions  of  the 
simulated  tropical  upper  tropospheric  troughs  in  our  model  agree 
quite  well  with  those  observed.   This  may  suggest  that  damping  due 
to  cumulus  transport  is  more  likely  to  account  for  the  vorticity 
sink  in  the  Tibetan  high  area.   V/e  also  found  singificant  transient 
wave  activities  at  the  upper  levels  in  both  experiments  and  they 
appear  to  be  drawing  energy  from  the  planetary  scale  barotropi- 
cally,  in  accordance  with  the  observation  by  Krishnamurti  (.1971). 
It  is  possible  that  they  play  an  important  role  in  reducing  the 
amplitudes  of  the  long  waves  which  agree  quite  well  with  those 
observed  despite  the  phase  problem.   This  point  is  borne  out  by 
comparison  with  the  excessively  large  amplitude  found  in  Holton 
and  Colton' s  linear  model  when  strong  damping  is  excluded.   In 
addition,  these  short  waves  are  most  active  in  the  fluctuating 
heating  case,  suggesting  that  local  barotropic  instability  at 


upper  levels  may  be  more  important  than  that  implied  by  the  time- 
mean  wind.   These  results  have  been  reported  by  Chang  and  Williams 
(1976b)  and  Chang  and  Pentimonti  (1977). 

Some  of  the  results  concerning  the  planetary-scale  monsoon, 
including  the  interpretation  of  the  stationary  equatorial  winter 
monsoon  circulation  as  a  viscous,  thermally-forced  Kelvin  wave 
(part  a) ,  the  preliminary  result  of  the  study  of  barotropic  in- 
stability of  the  upper- tropospheric  summer  monsoon  basic  flow 
(part  b) ,  and  the  preliminary  result  of  the  planetary  scale  vor- 
ticity  budget  using  forced  numerical  model  (part  c)  are  included 
in  a  review  paper  by  Chang  (1977b).   The  research  on  these  prob- 
lems is  being  continued  under  NSF  Grant  ATM77-14821  which  began 
on  15  September  1977. 
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ABSTRACT 

Solutions  to  the  wave-CISK  (conditional  instability  of  the  second  kind  with  cumulus  heating  being  in- 
duced by  low-level  internal  wave  convergence)  system  are  obtained  to  study  the  vertical  structure  of  margin- 
ally unstable  waves.  A  diabatic  heating  profile  is  specified  that  resembles  those  observed  and  those  theoret- 
ically derived  from  simple  parameterization  schemes.  Upper  and  lower  bounds  for  the  vertical  wavelength 
of  the  unstable  waves  under  normal  heating  conditions  are  established  through  analysis  of  the  frequency 
(stability)  equation.  The  lower  bound  excludes  the  possibility  of  e.xcitation  or  maintenance  of  short  vertical 
wavelengths  (relative  to  the  vertical  scale  of  heating)  by  wave-CISK.  The  calculated  growth  rates  indicate 
that  this  result  is  basically  insensitive  to  the  vertical  heating  profile.  The  vertical  structure  of  the  most 
unstable  waves  is  also  computed  and  the  possible  roles  played  by  CISK  in  large-scale  tropical  waves  are 
discussed  in  light  of  these  results. 


1.  Introduction 

The  possibility  that  large-scale  waves  with  short 
vertical  wavelengths  may  be  excited  in  the  tropical 
atmosphere  has  been  raised  by  Lindzen  (1967),  who 
showed  that  tfie  equatorially  trapped,  internal  gravity 
and  internal  Rossby  waves  usually  have  a  small 
positive  equivalent  depth  and  therefore  a  short  ver- 
tical wavelength.  Hoiton  (1969,  1972a)  later  pointed 
out  that  such  short  vertical  wavelengths,  if  they  exist 
in  the  troposphere,  would  cause  great  difficulties  for 
numerical  weather  prediction  in  the  equatorial  region 
because  of  the  stringent  requirement  for  high  vertical 
resolution.  Whether  short  vertical-scale  waves  actu- 
ally have  a  significant  presence  in  the  tropical  tropo- 
sphere or  not  is  unclear  due  to  inadequate  observations. 
Studies  based  on  conventional  radiosonde  data  (e.g., 
Wallace,  1971 ;  and  others)  tend  to  suggest  that  most 
of  the  variance  of  the  waves  in  the  tropical  tropo- 
sphere is  associated  with  vertical  scales  comparable 
to  the  scale  height.  However,  Madden  and  Zipser's 
(1970)  analysis  of  high-resolution  rawinsonde  data 
obtained  during  the  Line  Island  Experiment  showed 
that  a  quite  short  vertical  wavelength  (~3  km) 
existed  over  the  central  equatorial  Pacific.  Further 
observational  studies  thus  seem  necessary  to  resolve 
this  problem. 

Condensation  heating  due  to  deep  cumulus  con- 
vection has  been  recognized  as  a  primary  energy 
source  for  large-scale  tropical  motions.  Two  possible 
mechanisms  which  may  be  responsible  for  the  excita- 
tion of  large-scale  tropical  waves  by  heating  have 


been  investigated.  The  first  is  a  one-way  forcing  by 
the  cumulus  scale  which  receives  no  feedback  from 
the  large  scale.  Numerical  study  by  Hoiton  (1971, 
1972a,  1973)  and  analytical  work  by  Chang  (1976) 
have  shown  that  this  mechanism  excites  tropical 
waves  that  have  vertical  scales  comparable  to  those 
of  the  forcing,  which  are  usually  10-14  km.  The 
problem  of  short  vertical  wavelengths  thus  appears 
unimportant  for  forced  waves.  The  second  possible 
mechanism  is  the  conditional  instability  of  the  second 
kind  (CISK)  which  involves  a  feedback  to  cumulus 
convection  by  low-level  convergence  of  large-scale 
motion.  This  low-level  convergence  may  be  due  to 
frictionally  induced  cross-isobaric  flow,  as  in  the  case 
of  typhoons  (which  shall  be  called  "Ekman-CISK"), 
or  due  to  the  structure  of  the  internal  waves  (which 
shall  be  called  "wave-CISK").  The  Ekman-CISK  has 
been  shown  by  Chang  (1971)  and  Chang  and  Piwowar 
(1974)  to  possess  no  preferred  horizontal  scale  for 
synoptic  wave  disturbances.  Thus  Ekman-CISK  may 
not  explain  the  generation  of  synoptic-scale  waves, 
but  it  may  be  quite  important  in  maintaining  the 
waves  once  they  are  initiated.  The  wave-CISK  mecha- 
nism has  been  investigated  by  Hayashi  (1970),  Lindzen 
(1974)  and  Lindzen  et  al.  (1975).  Although  Lindzen 
(1974)  and  Lindzen  el  al.  (1975)  found  no  direct 
horizontal  scale  selection,  the  most  unstable  wave- 
CISK  mode  has  a  vertical  scale  seemingly  comparable 
to  the  depth  of  the  sub-cloud  layer  or  the  low-level 
convergent  moist  layer,  depending  on  the  particular 
model  aspects.  A  short  vertical  scale  may  therefore 
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be  important  for  the  excited  waves.  Using  the  cal- 
culated vertical  scale  they  also  speculated  that,  due 
to  the  Doppler-shifting  effect  of  the  variable  basic 
flow  and  the  small  growth  rates,  a  period  of  -^5  days 
associated  with  a  wavenumber  zero,  mixed  Rossby- 
gravity  mode  is  likely  to  dominate  the  tropical  spectra. 
This  periodicity  thus  provides  a  basic  frequency  for 
tropical  systems.  Depending  on  the  local  basic  flow, 
it  may  force  resonantly  higher  wavenumbers  that  may 
be  relevant  to  the  observed  waves. 

Since  the  growth  rates  associated  with  wave-CISK 
are  usually  very  small,  it  may  be  expected  that  the 
CISK  waves  have  properties  similar  to  neutrally 
forced  waves.  The  difference  in  the  selection  of  vertical 
scales  between  the  forced  models  and  the  wave-CISK 
models  thus  appears  to  be  worth  further  investigation. 
In  addition  to  its  importance  for  tropical  numerical 
weather  prediction,  the  vertical  scale  also  determines 
other  properties  of  the  excited  internal  waves.  The 
purpose  of  this  paper  is  to  examine  the  growth  rates 
associated  with  wave-CISK,  particularly  with  regard 
to  the  vertical  scale  selection  mechanism. 

2.  Basic  equations 

Neglecting  the  shear  of  the  mean  zonal  wind,  the 
linearized  zonal  and  meridional  momentum,  hydro- 
static, thermodynamic  energy  and  continuity  equa- 
tions for  a  single  zonal  wavenumber  k  on  an  equatorial 
beta-plane  may  be  written 


iu'U  —0yv  =  —  ik4), 

i(j:v+l3yu.= , 

dy 


d4> 
dz 


RT 
H 


iojT+ii'T  =  —, 

Cp 

dv  d 

iM-\ \-e''"—{e-'i"ib)  =0, 

dy  dz 


(1) 

(2) 

(3) 
(4) 
(5) 


where  ti,  v,  w,  T,  4>  and  Q  are  the  perturbation  zonal 
velocity,  meridional  velocity,  vertical  velocity,  tem- 
perature, geopotential  and  diabatic  heating  rate, 
respectively;  w  is  the  Doppler-shifted  frequency, 
H  a  constant  scale  height,  V  the  static  stability, 
Cp  the  specific  heat  at  constant  pressure,  R  the  gas 
constant,  /3  the  meridional  gradient  of  the  vertical 
component  of  earth's  vorticity,  y  the  meridional  co- 
ordinate, and  z=—H\n  (p/po)  is  the  vertical  coordi- 
nate with  p  the  pressure  and  po  a  reference  pressure. 
Eqs.  (l)-(5)  may  be  combined  into  a  single  equa- 
tion in  w  that  may  be  separated  into  meridional  and 


vertical  structure  equations  by  assuming  that 
ib  =  ZY..{y)wM)explz/{2H):\, 

n 

where  the  density  factor  exp  \^z/(2H)']  is  separated 
for  convenience.  The  meridional  structure  equation  is 


+  ( ;^2  +  _ 

df  \C0  ghn        ghn 


8Y\ 
F„  =  0, 


(6) 


where  g  is  the  gravitational  constant  and  hn,  the 
equivalent  depth,  is  the  separation  constant.  Matsuno 
(1966)  and  Lindzen  (1967)  have  shown  that  the  solu- 
tions to  (6)  which  satisfy  the  boundary  condition 


F„-^0    as     |y|->c 
and  lead  to  the  frequency  equation 


(7) 


\  u 


are 


igh^^ 


ghj     13 


=  2«-fl,     ^^=-1,0,  1,2,...,  (8) 
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1— (g/^.)' 


r     k 
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-Hn+li^) 


-ighnY^ 
O)  J 


Xexp(-^V2), 


where  ^=0Ks^'d~*y  and  the  Hermite  polynomials 
H„{^)  have  values  only  for  n'^0.  The  meridional 
velocity  solution  is  Vn(y)°^Hni^)  exp  (  —  ^^/l).  Lamb 
(1973)  has  solved  (6)  using  the  more  general  form 
of  the  confluent  hypergeometric  functions  without  the 
restriction  of  (7).  In  such  a  case,  the  lateral  boundary 
conditions  are  given  by  the  choice  of  vanishing  me- 
ridional velocity  at  some  finite  \y\  where  the  con- 
fluent hypergeometric  functions  have  a  node.  The 
Hermite  polynomials  are  a  sub-set  of  the  confluent 
hypergeometric  functions. 

The  vertical  structure  equation  which  is  most  rele- 
vant in  the  present  calculation  is 


dhVn  Qn 

dz^  ghn 


where 


\n'  = 

giln 
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1 

4H2 


(9) 


(10) 


(?;  = QnYrT'  exp[-2/(2i^)]. 
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The  parameter  Xn  is  a  measure  of  the  vertical  wave- 
number  and  is  generally  complex  for  unstable  waves. 
In  the  case  of  forced  modes  w  and  k  are  known  pa- 
rameters of  the  forcing  function  so  that  fu  and  X„  are 
obtained  from  (8)  and  (10)  respectively.  For  CISK 
modes,  on  the  other  hand,  X„  is  determined  as  an 
eigenvalue  of  (9)  and  w  is  then  found  from  (8)  and  (10). 
In  the  present  calculation,  the  heating  function  Q'  is 
assumed  to  possess  a  white-noise  distribution  for  all 
Hermite  modes  (all  values  of  n)  so  that  the  calcula- 
tion of  h  will  be  independent  of  n.  The  subscript  n 
will  therefore  be  dropped  in  the  subsequent  discus- 
sions. It  is  well  known  that  the  vertical  structure 
equation  (9)  can  also  be  derived  by  considering  two- 
dimensional  wave  motion  only. 

As  in  the  forced  model  by  Chang  (1976)   the  fol- 
lowing boundary  conditions  are  used  to  solve  (9) : 


w  = 


0, 


at    z  =  0 


where 


Cie'^'-\-Cie-'^',     at    2  =  Z( 


C  1  —  /"w  2 


(11a) 
(lib) 


Here  Zt  is  the  height  of  tropopause,  Ci,  Co  are  con- 
stants, and  the  condition  (lib)  results  from  the  re- 
quirement that  latent  heating  vanishes  at  Zt.  The 
parameter  r  is  a  reflection  coel'licient  which,  in  the 
absence  of  vertical  wind  shear,  is  given  by  the  speci- 
fication of  the  static  stability  distribution.  If  the 
heating  is  assumed  to  vanish  below  the  cloud  base  (sc), 
the  solution  to  (9)  may  be  written  in  the  form 


iX2_|_/,— 1X2 


U'{Z)=- 

zv{z)  =  - 


re'^'-'re  ' 
X(l+r) 
sinXz 


Q' 

sinXz — dz,     z^Zt, 


(12a) 


X(l-fr) 


{re'^'+e-'^')—dz 

gh  X(l+r) 

Q' 

X  /    sinXz—dz,     Zt>z>Zc,     (12b) 


gh 
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sInXz 
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{re^'+e-^'>'')—dz, 
gh 


The   solution   procedure   using   the   Green's   function 
technique  is  the  same  as  in  Lindzen  (1974). 

3.  The  heating  function 

The  parameterization  of  the  effect  of  condensation 
heating  due  to  cumulus  convection  is  a  major  problem 
in  tropical  modeling.  Important  contributions  on  this 
problem  have  been  made  by  the  sophisticated  schemes 
developed  by  Arakawa  and  Schubert  (1974),  Ooyama 
(1971),  and  others.  On  the  other  hand,  in  simplihed 


analytical  models,  simpler  schemes  must  be  used  to 
keep  the  mathematical  analysis  tractable.  Usually  in 
these  models  the  vertical  distribution  of  the  heating 
is  specified  as  a  given  function  that  is  invariant  in 
time.  The  major  criticism  of  this  type  of  parameteriza- 
tion is  that  the  results  are  sensitive  to  the  specified 
heating  profile  such  that  one  can  obtain  the  desired 
result  by  simply  tuning  the  heating  profile.  Another 
criticism  is  that  the  inherent  time-variation  of  the 
heating  function  due  to  the  complicated  scale  inter- 
action cannot  be  properly  included.  However,  recent 
observational  studies,  using  both  the  spectral  analysis 
technique  and  the  composite  technique,  have  produced 
a  consistent  picture  of  the  vertical  heating  profile  due 
to  latent  heat  release  over  the  tropical  oceanic  areas. 
The  profiles  observed  by  Nitta  (1970),  Wallace  (1971), 
Reed  and  Recker  (1971)  and  Williams  and  Gray  (1973) 
are  plotted  in  Fig.  1.  It  can  be  seen  that  in  each  case 
a  ma.ximum  occurs  in  the  middle  troposphere  between 
the  6  and  9  km  levels.  On  the  other  hand,  observa- 
tional analysis  of  the  divergence  and  vertical  velocity 
by  Schubert  and  Reed  (1975)  using  recent  GATE 
data  has  yielded  a  result  that  implies  a  vertical  heating 
profile  different  from  those  shown  in  Fig.  1.  They 
found  that  the  maximum  vertical  velocity,  which 
usually  coincides  with  the  maximum  heating,  based 
on  scaling  arguments  (Holton,  1972b),  occurs  at  a 
much  lower  level.  However,  analyses  of  the  same  data 
set  by  Petrossiants  ei  al.  (1975)  and  Antsipovich  et  al. 
(1975)  of  the  USSR  Hydrometeorological  Center  have 
resulted  in  profiles  which  are  in  close  agreement  with 
Fig.  1.  In  any  case,  the  profiles  included  in  Fig.  1  are 
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Fig.  1.  Vertical  heating  functions  deduced  from  observations  and 
from  simple  theoretical  parameterization  schemes. 
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based  on  a  much  larger  sample  of  data  and  may  be 
viewed  as  an  indication  that  the  tropical  atmosphere, 
on  a  statistical  basis,  tends  to  maintain  a  similar 
vertical  profile  for  the  large-scale  heating  function. 

Results  of  a  theoretical  calculation  by  Kuo  (1965) 
using  his  parameterization  scheme  are  also  included 
in  Fig.  1.  Kuo's  parameterization  was  originally  based 
on  the  assumption  of  horizontal  mixing  between  cloud 
air  and  environmental  air.  Recently  Kuo  (1974)  has 
shown  that  his  formulation  is  equivalent  to  the  con- 
sideration that  heating  is  due  to  the  adiabatic  com- 
pression of  the  descending  environment.  Kuo's  scheme 
has  been  used  in  man\-  tropical  numerical  models 
with  a  degree  of  success  not  inferior  to  any  other 
parameterization  scheme.  It  is  therefore  interesting  to 
compare  his  heating  profiles  with  those  observed.  In 
Fig.  1,  the  Kuo-1  curve  is  for  the  case  of  no  entrain- 
ment  and  the  Kuo-2  curve  is  for  the  case  of  moderate 
entrainment.  It  is  clear  that  the  latter  curve,  having 
a  maximum  near  7  km,  agrees  quite  well  with  those 
deduced  from  observations.  The  average  heating  profile 
obtained  from  the  GFDL  general  circulation  model 
as  calculated  by  Hayashi  (1973)  is  also  included, 
because  the  GFDL  model  uses  a  different  parame- 
terization scheme:  the  moist  convective  adjustment. 
Hayashi  decomposed  the  heating  function  into  com- 
ponents of  various  equatorial  beta-plane  modes.  Among 
them,  the  profile  for  the  Kelvin  waves  (;z=  — 1)  at 
the  equator  is  one  of  the  curves  least  resembling 
those  observed.  This  profile  is  also  plotted  in  Fig.  1. 
It  is  seen  that  even  for  this  highly  simplified  parame- 
terization scheme  the  maximum  heating  still  occurs 
near  6  km,  within  the  range  of  the  other  curves. 
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Q'(Z)  — 
Fig.  2.  The  specified  heating  function  Q'{z). 


Based  on  the  above  discussion  the  heating  functiorie 
for  the  wave-CISK  parameterization  is  specified  as 


Q'{z)  =  ^mNSiWhe"''  sin  -wz' ,     z^^z'^Zc 
(3'  =  0,  2>z,     or     z<Zc 


(13; 


where 


2   =- 


:-2, 


Az 


and  Az  —  Zi  —  Zc  is  a  measure  of  the  vertical  scale  ol 
the  heating.  The  coefficient  m  specifies  the  strength 
of  the  heating  and   the  factor  \   is  needed  because li' 
the  amplitude  of  the  Fourier  component  is  one-hal|E  ■ 
of   the   heating  maximum,   as  only  positive  conden 
sation  heating  is  permitted.  The  heating  is  proper-! 
tional  to  the  moisture  convergence  in  the  moist  layeief" 
which  is  represented  by  the  vertical  velocity  Wb  at  thelie 
top  of  the  moist  layer.  The  tropospheric  value  of  tWiay 
static  stability  Si  is  included  in  the  proportionalit} 
constant.  The  parameter  a  is  used  to  vary  the  maxi- 
mum heating  level  to  test  the  sensitivity  of  the  model 
The  coefficient  -V  is  a  normalization  factor  so  thatfett 
when  other  conditions  are  equal  the  total  amount  oiim 
heat  release  weighted  by  density  in  a  column  would' 
remain  the  same  with  different  values  of  a.  Thus, 


exp[-s/(2£?)]Js' 


e.xp[a2'— z/(2//)3  simrz'Jz' 


_2/i{l-exp[-Az/(2i/)]}{[a-Az/(2^)y+7r^} 

7r{exp[a-Az/(2//)]  +  l} 

Profiles  of   (13)   multiplied   by   the  density  factor 
exp  {z/2H)  are  plotted  in  Fig.  2  for  various  values" 
of  a.  The  range  of  —  l^a^O,  with  maximum  heating 
between  7  and  9  km,  appears  to  be  representative  of 
the  observed  and  theoretical  profiles  shown  in  Fig.  IJ 


4.  Analysis  of  the  stability  equation 

If  the  mixed  layer  with  top  at  Sb  (<Zc)  is  assume 
to   be    the   moisture   layer   which   provides   the   con 
vergence    for    CISK,    the   most    relevant   solution   in 
studying  the  instability  is  (12c)  evaluated  at  z^.  The 
solution  is 


:alt 


^Iv 


sinXzft  viWb  St  ttAz 


Wb=  — 


X      1+r  gh    2 


X 


\ 


where  m  =  mN  and  q^a—iXAz.  Equating  the  coefE' 
cients  of  wt  on  both  sides  and  using   (10)   leads  to 


n 


m 
E( 
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e  stability  equation 


9 


\        4//V 


[_{e''-\-l)c-'^''+r(e-''-\-].)e-'^"2 


X- 


(l+r)(^^=+7r-^) 


■1=0, 
(14) 


lere  we  have  assumed  tliat  (s'm  \z,,)/\~Zh.  It  can 
shown  that  q--\-!r'  =  0  is  not  a  solution  to  (14). 
J.  (14)  must  be  solved  for  the  complex  eigenvalue  X 
ing  either  numerical  or  graphical  methods.  However, 
fore  doing  so  some  valuable  information  regarding 
e  range  of  vertical  scale  of  the  unstable  solutions 
ly  be  obtained  by  analyzing  (14). 
The  complex  eigenvalue  X  may  be  written  as  X  =  Xr 
iX,-,  where  the  real  part  Xr  is  the  vertical  wave- 
imber  while  the  imaginary  part  X,  gives  the  growth 
te  of  the  waves.  If  the  positive  X^  is  chosen,  the  C2 
mponent  in  (lib)  represents  the  downward  phase 
upward  energy  propagation  and  the  Ci  component 
presents  the  upward  phase  or  downward  energy 
opagation.  It  then  follows  from  (10)  and  (8)  that 
X,  is  negative,  the  imaginary  part  of  the  frequency-  w 
11  be  negative  and  the  waves  will  grow  exponentially 
time.  Based  on  observations,  and  justified  a  pos- 
■iori,  we  shall  consider  only  the  case  of  marginally 
istable  waves,  or  very  small  growth  rates.  Thus  the 
nditions  applied  to  (14)  may  be  written 


X.2«X.2. 
addition,  we  shall  assume 


X.2»- 


1 
AW 


(15) 
(16) 


(17) 


lich  is  valid  for  all  practical  purposes  since  the 
ale  height  H  of  the  atmosphere  is  —7  km  and  (17) 
ily  places  a  restriction  of  ~88  km  for  the  maximum 
;rtical  wavelength. 

Eq.  (14)  will  now  be  analyzed  by  assuming  that 
=  0  and  that  the  stratosphere  {z>Zt)  and  troposphere 
^Z()  have  the  same  constant  static  stability  St. 
I  this  case  r  =  0  and  since  (f-\-iT-9^Q  the  stability 
[uation  may  be  re-written  as 

'(X)=>^(X)(r+7r2)=-w36Ac(x-+— Ve^+l) 
2  \        4//V 

Xexp(-fX,3,.)+(/--hx-  =  0. 

iparating  '^' {\)   into   its   real   and   imaginary   parts 


gives  two  equations: 

Re(^')=0, 
lm(^')=0. 

For  a  =  0,  (18)  may  be  expanded  to  give 
1 


(18) 


Re  (>!'')=    ^r--X 


4//2 


[e.xp(Xi3r)cosXrZ, 


+exp(X,2,.)cosX,s,]+2XrX.C^xp(X,z,)sinXrZ( 

2[7r--Az-'(X.2-Xr)] 

+exp(X.2,)sinX,s,]H =  0. 

TT/nz  bAz 

Simplifying    this   equation   with    the   conditions    (16) 
and  (17)  and  using  (15)  lead  to 

Re(^F') 

~  (X;-  cosXrZ,  — 2XrlX,!sinXrZ,)exp(—  |X,|2c) 

I 

+  (Xr-  cosXrZc  — 2X  J  X,  I  sinXrZ,.)exp(—  \  X,]  zj 


II 


2(7r'-Az2X.2) 


irrnzh^z 


=  0.      (19) 


III 


Since  '  cos  X^z,' ,  sinX^z,',  '  cos  X^z,. '  and  IsinX^z^! 
are  all  ^1  and  exp(— !  X,  j  z,.)<  1,  from  (16)  we  have, 
allowing  for  small  error  under  extreme  cases,' 


(20) 


|ii!^x,2| 

Thus  for  (19)  to  be  satisfied  we  may  write 

Max[Re(^')]  =  Max(I)+Max(II)+III>0,, 

where   Max  denotes   the  maximum   value   using   the 
range  (20) : 

Max(I)  =  Max(II)  =  XA 


Thus  we  have 


which  vields 


2X.2 


27r        2AzX.2 


mZiAz      7rmz6 

o 
TT" 

Xr-< , 

A2(Az— TrmZft) 
Az>xmz6. 


->o, 


(21) 


(22) 


'The  precise  ma\-imum  for  terms  I  and  II  is  |  j  (l-(-4e-)cos 
[tan~'(  — Se)  I  ]  e.xp  (  —  1 X,- 1  icJ !  Xr',  where  e=iX,:/Xri.  Inequality 
(20)  is  true  when  e  =  0.  For  e  =  0.1  which  is  our  assumption  the 
maximum  is  [1.02  exp(— IXj  jsjIXr^.  Obviously  if  !X,i  becomes 
large  the  e.xponential  factor  would  become  much  less  than  1. 
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Fig.  3.  Values  of  Re(^')  and  Inif*)  as  a  function  of  XrCkm"']  and  X,[km"'j.  The  crosses  in  the  Re(*) 
diagram  are  solutions  for  which  Re(*)  =  lm(*)  =  0. 


Choosing  the  minimum  values  (Min)  in  (20),  we  have 
Min[Re(^')]  =  Min(I)+Min(II)+III<0, 


so 


or 


giving 


2t        2A2\r2 

-2x.=+^: :^<o, 

mZbAz      irmzb 


Xr'>- 


Az{Az+TrmZb) 


(23) 


So  (21)  and  (23)  give  the  upper  and  lower  bounds 
for  \r.  It  is  interesting  to  note  that  Xr—^w/Az  in  the 
limit  of  m— >0,  corresponding  to  a  vertical  wavelength 
of  twice  the  heating  scale,  which  is  the  most  efficiently 
excited  wavelength  for  forced  waves  (Chang,  1976). 
However,  Xr  =  T/Az  cannot  be  a  solution  for  any 
finite  m  because  it  does  not  satisfy  (14). 

Now  we  shall  use  the  values  z,=  15  km,  2,.  =  1  km, 
Az=l4:  km,  26  =  0.4  km,  w  =  5.8  which  leads  to  m  =  9. 
The  value  of  m  is  calculated  by  fitting  (13)  to  Reed 
and  Recker's  (1971)  observed  heating  and  vertical 
velocity  profiles  for  the  wave  trough,  for  which  an 
estimate  of 

OT26==:2.32  km,  (24) 

can  be  made  assuming  that  6'(  =  1.22XlO~*  s~^  based 
on  a  tropospheric  value  of  r  =  3°C  km~^  Eq.  (24)  is 
a  good  approximation  as  long  as  S6<4  km,  so  the 
particular  choice  of  Zb  is  not  crucial.  From  Reed  and 
Recker's  data  Wb~0.29  cm  s~'  so  that  the  maximum 
heating  is  ~6.8°C  day"'.  Substituting  these  values 
into  (21)-(23)  we  obtain 


X  IT 

=  0.512>  |X,|  > =  0.167  km-i, 

6.14  18.8 


;,yj  the  vertical  wavelength,  L  =  2Tr/Xr,  is  limited  by 
the  following  bounds :  ' 

13  km<Z.<40km. 

Thus  waves  with  a  vertical  wavelength  shorter  than 
13  km  will  not  be  excited  by  wave-CISK.  If  such 
waves  are  first  excited  by  another  mechanism,  the 
above  result  indicates  that  they  cannot  be  maititained 
by  wave-CISK. 

5.  Eigenvalues  and  eigenfunctions 

a.  Growl  It  rales 


The  stability  equation  (14)  will  now  be  solved  using; 
the  graphical  method  with  the  accuracy  of  the  solu- 
tions improved  by  a  numerical  iteration  technique. 
The  graphical  method  is  illustrated  in  Fig.  3,  where 
the  real  and  imaginary  parts  of  ^(X)  given  by  (14) 
with  r  =  0,  a  =  0  and  m  —  6  are  plotted  as  a  function 
of  Xr  and  X,.  Here  the  constants  used  are  H  =  7  km, 
2,  =  15  km,  2,  =  1  km  and  26  =  0.4  km.  The  intersections 
of  the  zero  Re(^)  and  zero  Im(^)  lines  give  the 
eigenvalue  solution  X.  It  is  seen  that  there  exist  many 
solutions  with  the  most  (and  only)  unstable  solution 
having  the  smallest  X^  or  the  longest  wavelength.  The 
solutions  are  labeled  by  an  integer  k  with  increasing  k 
indicating  increased  X,  or  decreased  vertical  wavelength. 

The  four  longest  vertical  wavelength  solutions  are 
shown  in  Fig.  4  as  a  function  of  heating  strength  m. 
It  is  seen  that  |Xi|«X,  and  that  the  growth  rates  as 
represented  by  —X,  increase  with  m  while  the  vertical 
wavelengths  are  almost  constant.  This  diagram  also 
indicates  that  the  ^  =  1  solution  is  most  unstable 
throughout    the    range    m  =  4    to    w=10,    which    cor- 
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Fig.  4.  The  four  longest  vertical  wavelength  solutions  as  a 
function  of  the  heating  rate  m.  The  solid  curves  are  \,  and  the 
dash-dotted  curves  X;.  The  corresponding  vertical  wavelength 
Lt  =  2-K/\r  is  given  in  the  parentheses  of  the  Xr  coordinate.  The 
ratio  of  static  stabilities  between  stratosphere  and  troposphere 
is  1. 

responds  to  maximum  heating  rates  of  4.7  and  11.7°C 
day~^  if  ^6  =  0.29  cm  s~'.  In  fact,  it  is  the  only  un- 
stable solution  up  to  OT  =  9.  For  this  solution  the  value 
of  X  for  w  =  5.8  is  0.388-0.026^  [km"'],  which  is  well 
within  the  range  predicted  by  (21)  and  (23). 

The  sensitivity  of  the  eigenvalues  with  respect  to 
the  vertical  heating  prohle  (13)  is  shown  in  Fig.  5  in 
which  the  ^  =  1  to  ^  =  4  solutions  with  w  =  6  are  given 
as  a  function  of  a.  It  is  clear  that  the  vertical  wave- 
lengths are  all  insensitive  to  a.  For  the  growth  rates 
the  ^=1  solution  is  almost  independent  of  a,  but  the 
other  solutions  all  increase  with  a.  There  is  an  indica- 
tion that  as  a  increased  to  >5  the  shortest  vertical 
wavelength  becomes  the  most  unstable  solution.  How- 
ever, for  the  realistic  range  of  O^a^— 1,  which 
represents  the  statistical  mean  state  of  the  heating 
function  in  the  tropical  atmosphere,  ^=1  is  the  only 
unstable  solution. 

Solutions  for  a  different  static  stability  S^  in  the 
stratosphere  are  obtained  by  including  the  radiation 
condition  for  2>z,  as  the  third  boundary  condition. 
In  this  case  the  reflection  coefficient  is  derived  by 
continuity  requirements  of  w  and  ^  across  Zt : 

1-(V5,)^ 

r  = exp{-2iXtZc),     for    5,>5,. 

l  +  (5,/50^ 

The  result  for  Ss/St  =  2.5  is  shown  in  Fig.  6.   It  is 
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Fig.  5.  As  in  Fig.  4  except  that  the  solutions  are  given  as  a  function 
of  the  heating  profile  parameter  a. 

seen  that  the  general  conclusions  drawn  from  Fig.  5 
remain  valid,  although  waves  with  a  shorter  vertical 
wavelength  become  more  unstable  than  the  k—\  waves 
beginning  at  a  =  2  which  corresponds  to  a  maximum 
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Fig.  6.  As  in  Fig.  4  except  that  the  solutions  are  given  as  a 
function  of  the  heating  profile  parameter  a  and  that  the  ratio  of 
static  stabilities  between  stratosphere  and  troposphere  is  2.5. 
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heating  level  of  ~12  km  or  200  mb,  much  above  all     " 
observed  levels.  So  our  result  is  basically  insensitive 
to  the  vertical  heating  distribution. 

b.  Vertical  profiles 

The  vertical  structure  of  w,  corresponding  to  the 
X  for  k=  1  and  a=0  in  Fig.  6,  is  shown  in  Fig.  7  for  two 
heating  intensities  m=6  and  10.  It  can  be  seen  that 
the  tropospheric  amplitude  reaches  a  maximum  be- 
tween 6  and  9  km  and  has  a  minimum  at  the  tropo- 
pause  for  both  cases.  In  the  stratosphere  the  amplitude 
for  m  =  6  increases  at  a  slow  rate  with  height  due  to 
the  density  factor,  but  the  kinetic  energy  actually 
decreases  with  height.  This  decrease  is  due  to  the 
instability  of  the  waves  which  results  in  /\i<0.  The 
decrease  of  the  wave  kinetic  energy  is  much  faster 
for  the  more  unstable  w=10  case,  so  that  even  its 
amplitude  decreases  with  height.  The  vertical  phase 
distributions  for  both  cases  are  almost  the  same  and 
indicate  that  there  is  virtually  no  vertical  propagation 
in  the  troposphere  in  contrast  to  a  downward  propa- 
gation in  the  stratosphere.  The  vertical  wavelength 
computed  from  2ir/\r  is  16.2  km  in  the  troposphere 
and   10.2  km  in  the  stratosphere,  because  the  value 


*                                                \ 

W                                                               \ 
Ls=10.2km                         , 

Zt                                        \ 

v\ 

Lt  =16.2  km                                     \ 

nn=6   \ 

1 

\  m  =  IO 

/ 

Zc 

/ 

1          1          1 

1 

1  2  3 

AMPLITUDE 


O 
PHASE 


Fig.  7.  Amplitude  and  phase  of  w  as  a  function  of  height  for 
;«  =  6  (solid  lines)  and  w=10  (dotted  lines).  The  phase  lines  for 
both  cases  are  approximately  the  same. 
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Fig.  8.  As  in  Fig.  7  except  for  T. 


of  Xr  is  increased  in  the  stratosphere  by  approximately 
a  factor  of  (Ss/Si)'',  as  can  be  seen  from  (10). 

Vertical  structures  of  other  variables  can  be  com- 
puted utilizing  the  basic  equations  (l)-(5)  and  the 
solutions  to  the  meridional  structure  equation  (6). 
Here  only  the  temperature  structure  for  the  Kelvin 
wave  case  (»=  — 1)  will  be  shown.  From  (4)  and  (8) 
we  have 

RQ 


r„=_i=  — 


ilwS )// 

\  H  cj 

kR{gh)i 


(25) 


It  is  obvious  from  (25)  that  the  amplitude  of  the 
temperature  perturbation  decreases  as  the  zonal  wave- 
number  k  increases.  The  Tn=-i  structure  corresponding 
to  the  w  in  Fig.  7,  is  shown  in  Fig.  8.  It  can  be  seen 
that  the  amplitude  distribution  is  quite  sensitive  to 
the  strength  of  the  heating.  The  discontinuity  at  2,  is 
due  to  the  layered  static  stability  distribution  assumed 
in  our  model.  The  large  temperature  perturbations  in 
the  stratosphere  are,  of  course,  due  to  the  large  static 
stability  there.  Figs.  7  and  8  are  constructed  using 
nondimensional  units.  The  two  figures  may  be  com- 
pared  by   letting   the   unit   for   the  w  amplitude   be 
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cm  s~' ;  in  this  case  the  corresponding  unit  for  the 
amplitude  will  be  a~'  C°C3,  where  a  =  kX  (radius  of 
he  earth)  is  the  number  of  waves  around  the  equator. 
"or  these  units  the  amplitude  of  w  is  quite  close  to 
he  typically  observed  values  of  "matured"  tropical 
vaves,  with  maximum  ^3  cm  s~^  for  m  =  6  and 
4.2  cm  s~'  for  m  —  10.  On  the  other  hand,  the  maxi- 
num  T  perturbation  in  the  troposphere  is  almost 
LO°C  for  wavenumber  (a)  1  and  even  larger  in  the 
stratosphere.  Lindzen  (personal  communication)  has 
suggested  that  these  enormously  large  temperature 
fluctuations  indicate  that  large  damping  effect  of 
vertical  momentum  transport  due  to  cumulus  con- 
vection must  be  important.  However,  for  shorter  waves 
the  temperature  perturbation  is  proportionally  smaller. 
If  we  consider  a  zonal  wavelength  of  4000  km,  then 
a~10  and  the  temperature  fluctuations  are  all  <1°C 
in  the  troposphere,  in  good  agreement  with  observa- 
tions and  with  Holton  (1971)'s  numerical  calculations. 
In  fact,  this  magnitude  estimate  may  explain  why 
atmospheric  Kelvin  waves,  which  have  a  T  amplitude 
~2°C  in  the  lower  stratosphere,  are  not  observed  in 
the  troposphere.  For  such  waves  the  vertical  velocity 
in  the  troposphere  would  be  an  order  of  magnitude 
smaller  than  for  the  easterly  waves  which  have  a 
typical  zonal  wavelength  ~3000  km  and  a  maximum 
vertical  velocity  ^3  cm  s~'.  On  the  other  hand,  damp- 
ing effect  of  a  certain  magnitude  must  pla\-  a  role  for 
the  easterly  waves  to  explain  the  complete  absence 
of  them  in  the  stratosphere.  Although  these  estimates 
are  based  on  Kelvin  waves  only,  for  other  wave  modes 
a  relationship  between  the  tropospheric  temperature 
amplitude  and  the  zonal  scale  must  also  exist  at 
latitudes  somewhat  away  from  the  equator  because 
of  the  following  scaling  arguments.  Wallace  (1971) 
and  Holton  (1972b)  have  shown  that,  for  weak, 
tropical,  synoptic-scale  motions  with  a  depth  scale 
comparable  to  the  scale  height,  a  diabatic  heating 
rate  of  5-10°C  day~^  is  almost  completely  balanced 
by  the  adiabatic  cooling.  Thus  the  temperature  fluc- 
tuation is  one  order  of  magnitude  smaller  than  either 
term.  This  balance  may  be  attributed  to  the  smallness 
of  the  Rossby  number  which  is  ~0(1)  or  larger  for 
tropical  synoptic-scale  motions.  For  planetary-scale 
motions  the  length  scale  is  one  order  larger  so  that 
the  Rossby  number  is  -^0(10"')  at  intermediate 
tropical  latitudes,  and  the  scale  analysis  indicates  that 
the  pressure  and  temperature  perturbations  must  be 
increased  by  one  order  of  magnitude. 

The  phase  diagrams  in  Figs.  7  and  8  indicate  that 
in  the  lower  troposphere  and  the  entire  stratosphere 
the  temperature  leads  the  vertical  velocity  by  ~^  cycle. 
In  the  middle  and  upper  troposphere  (7-12  km  for 
w  =  6  and  5-12  km  for  ot=10),  they  are  in  phase 
except  near  the  tropopause  (12-14  km),  where  they 
are  out  of  phase.  These  phase  relationships  are  con- 
sistent with  the  energetics  because  the  conditions  for 


wave  growth  are 


Jo 


w'T'dz>0, 


Q'T'dz>0. 


From  (25)  and  (10)  it  can  be  shown  that  these  con- 
ditions are  satisfied  when  X,<0  and 


SRe{w)dz<R       Qdz/(Hcp) 


The  latter  condition  implies  that  the  diabatic  heating 
rate  exceeds  the  adiabatic  cooling  rate  which,  of 
course,  is  responsible  for  the  continuous  growth  of  the 
waves  with  time. 

c.  Equatorial  (rapping  scales 

The  meridional  scale  of  the  unstable  waves  may  be 
estimated  from  the  Hermite  solutions  of  the  meridional 
equation  (6).  The  north-south,  e-folding  width  is 

y.  =  Re[2(g^)V/3]s 

and  the  maximum  amplitude  for  the  mixed  Rossby- 
gravity  mode  {n  =  0)  occurs  at 

The  maximum  for  the  Kelvin  mode  occurs  at  the 
equator.  From  results  obtained  for  the  k  =  l  waves, 
\~O.388-0.026z  [km-']  so  that  (ghY'-lB.S-lJi  Im 
s-']  from  (10).  It  follows  that  ^^-1480  km  and  yo 
-1050  km  with  /i  =  2. 21X10-"  m"'  s"'.  Both  of  these 
are  reasonable  values  in  terms  of  the  trapping  of 
waves  on  the  equatorial  beta-plane. 

6.  Concluding  remarks 

We  have  shown  that  tropical  waves  cannot  remain 
unstable  with  respect  to  wave-CISK  if  the  vertical 
wavelength  is  much  smaller  than  the  vertical  scale 
of  heating.  This  result  holds  for  a  number  of  reason- 
able heating  profiles  that  have  been  observed  or 
theoreticall}'  derived.  The  dispersive  relationship  of 
the  growth  rates  of  various  equatorial  beta-plane 
modes  has  been  shown  by  Lindzen  (1974)  to  be  un- 
satisfactory in  selecting  the  horizontal  scale  of  tropical 
waves.  From  the  computed  eigenvalues  it  can  be 
shown  that,  even  under  the  extremely  intense  heating 
case  of  w=10,  the  growth  rates  for  Rossby  waves, 
mixed  Rossby-gravity  waves,  and  the  longer  zonal- 
scale  Kelvin  waves  are  all  very  small.  This  result  is 
similar  to  that  obtained  by  Chang  and  Piwowar 
(1974)  for  the  Ekman-CISK  and  indicates  that  wave- 
CISK  is  not  responsible  for  the  initial  selection  mecha- 
nism of  the  tropical  waves.  Once  the  waves  are  e.xcited 
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by  some  other  mechanism,  however,  the  internal  wave 
convergence  could  induce  organized  cumulus  convec- 
tion and  release  wave-CISK  as  an  energy  source.  The 
present  calculation  indicates  that  under  such  cases 
the  shorter  vertical-scale  waves  will  still  not  be  main- 
tained by  wave-CISK.  There  is  no  definite  explana- 
tion to  the  fact  that  the  gravity  modes  and  the  shorter 
zonal  scale  Kelvin  modes,  which  are  usually  not 
observed  in  the  tropical  atmosphere,  are  most  unstable. 
However,  these  waves  have  very  short  periods  and 
various  arguments  have  been  given  against  their  oc- 
currence on  a  regular  basis.  For  example,  if  one  con- 
siders the  lifting  process^  of  CISK  (Ooyama,  1969; 
Holton,  1972b),  the  short  periods  will  be  insufficient 
for  a  low-level  air  parcel  to  reach  its  condensation 
level.  In  any  case,  the  result  that  wave-CISK  can 
neither  excite  nor  maintain  waves  with  short  vertical 
wavelengths  remains  valid  independent  of  these  con- 
siderations. The  concern  about  short  vertical-scale 
excitation  due  to  wave-CISK  in  a  tropical  numerical 
model  may  thus  be  relieved.  The  resonant  forcing  of 
tropical  waves  by  a  fundamental  wave-CISK  fre- 
quency associated  with  the  zonally-symmetric,  mixed 
Rossby-gravity  mode  as  hypothesized  by  Lindzen 
(1974)  also  does  not  appear  to  be  plausible,  because 
the  equivalent  depth  would  be  too  large  to  give  a 
fundamental  periodicity  near  5  days. 

The  calculated  vertical  structure  of  w  and  T  applies 
to  all  waves  maintained  by  heating  and  is  not  re- 
stricted to  wave-CISK.  The  eigenvalue  X  found  by 
the  wave-CISK  equation  used  here  may  alternatively 
be  viewed  as  being  excited  by  a  specified  forcing 
[as  in  the  model  by  Chang  (1976)]  with  a  small 
linear  damping.  The  damping  coefficient  would  then 
appear  in  the  complex  Doppler-shifted  frequency  in 
the  same  way  an  an  imaginary  frequency  due  to 
instability.  The  computed  magnitudes  of  w  and  T 
suggest  that  the  observed  planetary-scale  Kelvin 
waves  are  likely  to  be  a  result  of  forcing  rather  than 
instability,  because  a  small  forcing  function  can  ex- 
plain the  amplitudes  in  the  stratosphere  and  tropo- 
sphere. However,  an  instability  consistent  with  the 
present  CISK  parameterization  should  eventually 
amplify  the  waves  to  a  finite  amplitude  w  comparable 
to  that  observed  for  large-scale,  cumulus-associated 
tropical  motions,  but  with  a  temperature  amplitude 
too  large  to  be  realistic.  The  small  low-level  vertical 
velocities  of  these  waves  is  also  an  indication  that 
CISK  is  unlikely  to  be  important.  On  the  other  hand, 
the  tropospheric  amplitude  of  the  synoptic-scale  waves 
may  quite  possibly  be  supported  by  a  CISK-type 
process.  These  waves  have  a  wavelength  ~  2000-4000 
km  and  a  period  ~5  days  which  lead  to  a  phase  speed 
relative  to  the  ground  of  ~5-9  m  s~'.  Such  a  phase 


-Lindzen  el  al.  (1975)  have  alternatively  considered  that  the 
CISK  is  not  a  lifting  process  but  a  low-level  organizing  process 
which  does  not  require  lifting. 


speed  would  make  the  waves  greatly  attenuated  before 
reaching  the  stratosphere,  because  the  Doppler-shifted 
phase  speed  is  very  small  near  and  below  the  tropo- 
pause  and  the  damping  can  become  very  efficient 
(Lindzen,  1971). 
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ABSTRACT 

The  problem  of  scale-selection  of  Kelvin  waves  in  the  stratosphere  by  forcing  from  tropospheric  heating 
is  analyzed  using  a  simple  linear  model.  The  effect  of  vertical  wind  shear  is  excluded  because  the  phase 
speed  of  the  waves  is  fast  relative  to  the  range  of  the  mean  zonal  wind  in  the  vicinity  of  the  tropopause  at 
which  level  the  upward  energy  flux  due  to  forcing  is  evaluated.  Results  of  this  analysis  modify  Hoi  ton's 
( 1973)  theory  in  that  1 )  the  forcing  is  most  efficient  for  the  longest  zonal  wavelength  even  if  the  heat  sources 
are  distributed  randomly,  and  2)  the  most  favored  vertical  wavelength  of  the  excited  waves  is  about  twice 
the  vertical  scale  of  heating.  The  calculated  vertical  wavelengths  exceed  slightly  those  observed  and  the 
discrepancies  are  discussed. 


1.  Introduction 

There  is  ample  evidence  on  the  existence  of  atmo- 
spheric Kelvin  waves  in  the  equatorial  stratosphere 
(Wallace  and  Kousky,  1968;  and  others).  These  waves 
propagate  eastward  and  downward  with  a  zonal  wave- 
number  of  1-2,  a  periodicity  of  10-15  days,  and  a 
vertical  wavelength  of  8-12  km.  Their  zonal  wind 
component  has  a  typical  amplitude  of  10  m  s~"'  at  the 
equator  while  no  appreciable  fluctuation  of  the  me- 
ridional wind  component  is  observed.  The  zonal  wind 
and  pressure  perturbations  are  both  symmetric  with 
respect  to  the  equator  and  are  close  to  geostrophic 
balance  in  the  meridional  direction,  but  the  waves 
behave  like  internal  gravity  waves  in  the  zonal  and 
vertical  directions.  More  detailed  observational  and 
theoretical  descriptions  of  these  waves  are  given  in 
Wallace  (1973)  and  Holton  (1975).  These  waves  are 
believed  to  play  a  very  important  role  in  the  quasi- 
biennial  oscillation  of  the  equatorial  stratosphere  by 
supplying  westerly  momentum  to  the  zonal  mean  flow 
(Lind'zen  and  Holton,  1968;  Holton  and  Lindzen, 
1972). 

The  vertical  structure  of  the  Kelvin  waves  indicates 
that  they  are  maintained  by  energy  propagating  from 
the  troposphere.  In  addition,  the  possibility  that  the\- 
are  forced  by  interaction  with  mid-latitude  motions, 
as  suggested  by  Mak  (1969),  is  quite  small  because 
of  the  small  meridional  wind  perturbations  which  rule 
out  the  importance  of  lateral  energy  fluxes.  Dynamic 
instability  also  appears  unlikely  due  to  the  very  long 
zonal  scale  and  the  fast  westerly  phase  speed  of  the 
waves. 

Due  to  the  presence  of  the  large  amount  of  deep 
cumulus  convection  in  the  tropical  troposphere,  it  has 


been  suggested  that  a  plausible  energ\-  source  fo 
these  waves  is  the  latent  heat  released  by  cumulu 
towers.  There  are  two  possible  mechanisms  by  whic 
the  latent  heating  in  the  troposphere  may  maintai 
the  Kelvin  waves : 

1)  A  two-way  interaction  such  that  the  cumuli 
convection,  which  supplies  energy  to  the  waves, 
itself  controlled, by  the  large-scale  wave  motion  fiek 
This  scale  interaction  results  in  an  unstable  situatic 
which  is  generally  called  the  conditional  instability  > 
the  second  kind  (CISK). 

2)  The  waves  are  forced  by  cumulus  heating  whic; 
is  controlled  by  processes  unrelated  to  the  waves 
Havashi  (1970)  and  Lindzen  (1974)  have  investigate! 
the  first  mechanism  and  found  that  the  growth  ratt 
of  CISK  are  unsatisfactory  in  explaining  the  Kelvi 
waves.  Holton  (1972,  1973)  studied  the  second  meek 
nism  with  a  numerical  diagnostic  model  in  which  th 
diabatic  heating  due  to  cumulus  convection  is  spec 
hed.  Damping  in  the  form  of  Rayleigh  friction  an 
Newtonian  cooling  are  also  included.  He  found  thj 
the  wave  response  in  the  stratosphere  to  a  trope 
spheric  heat  source  resembling  the  Kelvin  wave  moc 
has  a  structure  in  close  agreement  with  observation: 
However,  the  observed  periodicity  of  the  Kelyi 
waves  must  be  specified  although  no  correspondin 
peak  in  the  cloud  brightness  spectra  can  be  foum 
On  the  other  hand,  his  numerical  calculations  suggei 
that  the  vertical  scale  of  the  e.xcited  waves  is  clos 
to  that  of  the  specified  forcing.  In  equatorial  wav 
theory  (Holton  and  Lindzen,  1968)  the  vertical  sea 
determines  the  wavenumber-frequency  relationshij 
Bv  postulating  a  red-noise  spectral  distribution  of  tW 
tropospheric  heat  sources,  Holton  suggested  that  th 
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vertical  scale  and  the  large  low-frequency  variance 
of  the  forcing  result  in  a  band-pass  selectivity  of  the 
Kelvin  waves. 

Helton's  theory  may  indeed  be  valid,  especially 
because  an  oscillation  of  the  monsoon  heating  in  South 
Asia,  with  a  quasi-periodicity  of  about  15  da>s,  has 
been  observed  during  the  northern  summer,  although 
it  is  not  clear  whether  this  periodicity  exists  in  other 
seasons.  The  purpose  of  this  paper  is  to  point  out 
that  the  atmosphere  acts  like  a  band-pass  filter  even 
if  the  heat  sources  of  the  troposphere  are  distributed 
randomly  in  the  frequency  domain.  The  relationship 
between  the  vertical  scale  of  the  response  and  that 
of  heating  will  also  be  examined. 

2.  The  forced  model 

The  zonal  mean  wind  in  the  tropical  stratosphere 
has  considerable  vertical  shear  which  may  influence 
the  upward  propagation  of  wave  energy  (Lindzen, 
1971).  However,  due  to  the  fast  phase  speed  of  the 
Kelvin  waves  the  shear  is  probably  not  very  signifi- 
cant at  the  tropopause  level  for  which  our  discussion 
will  be  most  relevant.  The  zonal  mean  wind  will 
therefore  be  neglected  in  our  simple  analysis.  The 
effect  of  damping  will  also  be  excluded.  The  linearized 
zonal  momentum,  meridional  momentum,  hydrostatic, 
thermodynamic  energy  and  continuity  equations  on 
an  equatorial  beta-plane,  with  a  single  zonal  wave- 
number  k,  a  Doppler-shifted  phase  speed  c,  and  no 
perturbation  in  the  meridional  velocity  may  then  be 
written 

(1) 


(2) 
(3) 
(4) 
(5) 


Kelvin  waves : 


ikcH  = 

■■-ik4>, 

l5yu  = 

d<t> 
dy 

dz 

RT 

h' 

Q 

-ikcT-\-wT=—, 


iku+e''"—{e-'i"w)=0, 
dz 


where  u,  w,  T,  4>  and  Q  are  the  perturbation  zonal 
velocity,  vertical  velocity,  temperature,  geopotential 
and  the  diabatic  heating  rate,  respectively;  //  is  a 
constant  scale  height,  T  the  static  stability,  Cp  the 
specific  heat  at  constant  pressure,  R  the  gas  constant, 
iS  the  meridional  gradient  of  the  vertical  component 
of  earth's  vorticity,  y  the  meridional  coordinate,  and 
z=  —H  \n(p' po)  is  the  vertical  coordinate  with  p  the 
pressure  and  po  a  reference  pressure. 

Eqs.    (l)-(2)   give  the  meridional  structure  of  the 


I3y' 
w,  (/>ocexp( ), 


(6) 


where  the  condition  c>0  is  required  to  satisfy  the 
trapping  condition  on  an  equatorial  /S-plane.  If  a 
height-dependence  factor,  exp{z/2H),  is  separated  from 
the  perturbation  quantities,  Eqs.  (1)  and  (3)-(5)  may 
be  combined  to  form  a  single  equation  in  w: 


d-w' 


-xv=- 


Q' 


dz^ 


(7) 


where 
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S 
R 


Q'- 


CpH 


_Qg-./(2//) 


(8) 


R 

5=-r 

H 


The  parameter  X  is  a  measure  of  the  vertical  wave- 
number.  If  the  heating  function  Q'  is  given,  (7)  can 
be  solved  with  suitable  boundary  conditions.  In  our 
problem  the  following  boundary  conditions  are  used : 


where 


0, 

at     s  =  0 

(9a) 

Cie^^'+C^e-'^^ 

at     z  =  Zi 

(9b) 

C\  =  rC2. 


Here  z,  is  the  height  of  tropopause  and  the  condi- 
tion (9b)  results  from  the  requirement  that  latent 
heating  vanishes  at  z^.  The  parameter  r  is  a  reflection 
coefficient  which,  in  the  absence  of  vertical  wind  shear, 
is  given  by  the  model  specification  of  the  static  sta- 
bility distribution.  The  solution  to  (7)  at  z  =  Zt  may 
now  be  written  as 


w{zt)  =  — 


re 


i\zt 


+e- 


-i\zi 


X(l+r) 
where  Zc  is  the  height  of  the  cloud  base 

3.  The  energy  flux  at  tropopause 


Q' 

sin  {\z)—dz 


(10) 


We  will  consider  the  vertical  energy  flux  4>'w\_ 
where  the  overbar  denotes  the  zonal  average  and 
<j>'  =  <pexp{  —  z/2H),  at  the  tropopause  level  (z,),  to  be 
a  measure  of  the  efficiency  of  forcing.  Despite  the 
actual  existence  of  the  vertical  shear  of  the  zonal 
mean  flow  in  the  stratosphere,  the  response  of  the 
waves  should  be  largelv'  determined  by  this  quantity. 
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At  z,  the  solution  is  given  by  (9b)  and  only  the  C2 
term  needs  to  be  considered  for  upward  energy  propa- 
gation or  downward  phase  propagation. 
Eqs.  (1)  and  (5)  may  now  be  rewritten  as 

—  ikcu'=  —ik(j>', 


iku'-{i\-\ V  =  0, 

\        2HJ 


where  z«'  =  m  exp(  — 2/2/7).  From  these  relationships  it 
follows  immediately  that 


(j,'  =  cu'  =  -(x W. 


Since  E~l  km  for  the  atmosphere,  X»l/(2//)  if  the 
vertical  wavelength  «88  km.  This  is  the  case  for  all 
practical  purposes  so  we  may  write 

c\ 

(t>'~—w'. 
k 

From  the  definition  of  \  given  in  Eq.  (8)  we  may  also 
write 


which  leads  to 


and 


k 

5i_ 

<P  "W   ^ W  '. 

k 


(11) 


Eq.  (11)  indicates  that,  if  the  amplitude  of  the  ver- 
tical velocity  perturbation  at  the  tropopause  level  is 
fixed,  the  upward  energy  flux  will  be  inversely  propor- 
tional to  the  zonal  wavenumber.  Examination  of  (7) 
reveals  that  w'  depends  only  on  X  (c  is  uniquely  de- 
termined by  X)  as  long  as  Q  is  independent  of  k.  Thus 
for  given  vertical  scale  the  forcing  at  the  tropopause 
level  becomes  maximum  for  the  longest  horizontal 
wavelength.  This  result  obviously  remains  the  same 
if  one  considers  the  upward  momentum  flux  ti'w', 
because  ii'  =  <i>'/c. 

In  order  to  examine  the  selection  of  vertical  scales, 
we  assume  a  heating  function 


mSq(y)  sirnri 

lo. 


z>Zt  or  z<Zc 


(12) 


where  m  specifies  the  strength  of  the  heating  and 
Az  =  Zt  —  Zc  is  a  measure  of  the  vertical  scale  of  the 
heating.  The  function  q(y)  which  specifies  the  me- 
ridional  distribution  of  heating   will   be   assumed    to 


take  the  form  of  (6)  for  the  time  being.  The  vertic 
profile  as  given  by  (12)  resembles  closely  that  of  tl 
typical  observations  (Chang,  1976).  For  the  simple 
case,  we  shall  first  assume  that  the  static  stabilil 
parameter  5  has  a  constant  tropospheric  value 

throughout  the  atmosphere.  In  this  case  r  =  0  and  onj 
upward  energy  propagation  at  z  =  Zi  is  allowed.  Tlj 
solution  (10)  at  Zt  is  now  found  to  be 


w'  =-w  e.xp[  — zX2(  — /3_yV(2c)],     2  =  z<, 


where 


W7rAsX2+ 1/(4^)2       rX 
w= sin  -{zi+z,) 


X 


■X-Az- 


L2 


-s(y) 


(i; 


(M 


denotes  the  amplitude  of  the  solution  in  this  simp 
model.  We  note  here  that  no  actual  singularity  exis 
in  (14).  The  vertical  energy  flux  at  the  equator  anj 
z  =  2(  is  thus 

5S_ 

<t>'w'  =  —w-. 
2k 

From  (14)  if  one  assumes  X»l/(4//-),  it  can  be  show 
that  this  energy  flux  reaches  maximum  at 


i 


X  =  - 


Az 


where 


ib^^ sin- 

16 

"    Tf                        "1 

(2,+Zc)     . 

_2Az             J 

The  tropospheric  forcing  is  therefore  most  efficient  a 

a  vertical  wavelength 

27r 
L.s  — =  2Az, 
X 

or  twice  the  vertical  scale  of  the  heating.  This  resul' 
is  somewhat  similar  to  those  of  Green  (1965)  anc 
Lindzen  (1966).  A  plot  of  the  energy  flux  profile  at  Zj 
is  given  in  Fig.  1  in  which  it  is  assumed  that  Z(  =  15j 
and  H  =  7Zc.  If  we  measure  the  vertical  scale  of  Hoi 
ton's  (1972,  1973)  heating  function  the  value  Az  car 
be  found  to  be  ~11  km,  which  predicts  a  maximun 
response  of  22  km  vertical  wavelength  in  the  tropoJ 
sphere.  Since  in  reality  the  stratospheric  static  sta- 
bility is  about  three  times  that  of  the  troposphere- 
using  (8)  we  may  estimate  that  the  stratospheri(i 
vertical  wavelength  corresponding  to  a  22  km  tropo- 
spheric value  would  be 

Z,«Li  — j     =22X3-'^12.8km, 

where  the  subscript  .J  denotes  stratosphere.  If  wf 
assume  Az—  14  km,  then  L^—  16.3  km. 
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It  may  be  of  some  interest  to  consider  the  effect 
of  different  static  stability  between  stratosphere  and 
troposphere  on  the  solution  (13),  because  the  "effi- 
ciency" of  reflection  of  wave  energy  at  z  =  Zt  may  be 
scale-dependent.  In  this  case  we  have 

S  =  S„     for  z>Zt, 
S=Si,     for  z^Zi. 

Another  boundary  condition  is  now  needed  and  the 
simplest  form  available  is  the  radiation  condition 
above  z,,  i.e., 


w 


■■A  exp[-iX,2-/33;V(2c)],     for  z>2,,      (15) 


where  X,  is  the  vertical  wavenumber  in  the  strato- 
sphere. The  continuity  conditions  of  w'  and  </>'  across 
z,  lead  to 


exp{  —  2iXtZi),     for     S,>St, 


where  X(  is  the  vertical  wavenumber  in  the  tropo- 
sphere. The  solution  (10)  at  z,  in  the  form  of  (15)  is 


.4=- 


2Xt  exp[/(X,  — Xt)Z(]tv 
(X.+X,)-(X,-X,)  exp(-2/X,z,)' 


where  w  is  given  by  (14)  with  ail  X  replaced  by  X,. 
The  vertical  energy  flux  at  the  equator  and  z,  is  now 

2k  2k 

where 

£  =  2[(l+5,/5,)-f(l-5./5,)cos2X(Z,]-' 

is  the  effect  of  different  static  stabilities  above  and  below 
Z(  and  the  asterisk  denotes  complex  conjugate.  For 
Ss=3St,  E  varies  between  a  maximum  of  1  when 
X«=0,  Tr/zt,  2ir/zt,  . .  .,  and  a  minimum  of  0.33  when 
Xt=7r/(2z«),  37r/(2z(),  5Tr/(2s(),  ....  Assuming  Zj=  15 
km  and  Az=  14  km,  the  variation  of  E  is  plotted  in 
the  upper  section  of  Fig.  1  and  the  modified  energy 
flux  profile  is  indicated  by  the  dashed  curve  in  the 
lower  section.  It  is  seen  that  the  most  efficient  forcing 
is  slightly  larger  than  twice  the  scale  of  heating.  Hence 
we  may  conclude  that  the  variation  in  static  stability 
between  troposphere  and  stratosphere  does  not  cause 
great  effect  in  the  vertical  scale  of  the  most  favored 
response. 

4.  Concluding  remarks 

We    have    shown    that    for    randomly    distributed 
tropospheric  heat  sources   the   Kelvin  waves   of   the 


0W 


4H, 


Fig.  1.  Profile  of  vertical  wave  energy  flux  <t>'w'  as  a  function  of 
vertical  wavelength  Lz-  The  solid  curve  in  the  lower  diagram  is 
for  the  case  S,  =  St  and  the  dashed  curve  for  the  case  S,  =  3Sf  The 
upper  diagram  is  the  effect  on  E  of  the  difference  in  static  stability 
on  the  vertical  wave  energy  flux  profile  for  the  case  S,  =  3Si. 


largest  zonal  scale  are  e.xcited  most  efficiently  in  the 
absence  of  damping.  This  may  e.xplain  why  most  of 
the  Kelvin  waves  observed  in  the  lower  stratosphere 
are  of  zonal  wavenumber  1  and  some  are  of  wave- 
number  2.  The  vertical  wind  shear  has  been  excluded 
from  our  analysis  but  it  should  not  significantly  alter 
this  conclusion. 

The  finding  that  the  most  favored  vertical  wave- 
length is  about  twice  the  vertical  scale  of  heatmg 
gives  a  reasonable  estimate  for  Kelvin  waves  in  the 
lower  stratosphere,  although  the  calculated  wave- 
length of  13-16  km  is  slightly  larger  than  those  ob- 
served. Several  influences  which  are  excluded  in  our 
model  may  alter  this  aspect.  For  example,  Lindzen 
(1971)  found  that  the  focusing  effect  of  the  vertical 
wind  shear  in  the  stratosphere  would  reduce  the 
vertical  scale  of  the  waves.  A  two-scale  expansion 
technique  can  be  similarly  applied  to  our  forced  model 
and  the  resultant  prevailing  vertical  wavelength  in 
the  stratosphere  should  be  reduced. 

Another  effect  which  has  been  excluded  is  the 
actual  north-south  distribution  of  the  heating  func- 
tion. The  meridional  structure  equation  (6)  indicates 
that  the  meridional  scale  of  the  waves  is  coupled  with 
the  vertical  scale  through  the  phase  speed  c.  The 
variation  of  the  meridional  scale  with  the  vertical 
scale  may  be  considered  slow  because  the  e-folding 
width  of  the  Gaussian  distribution  is  proportional  only 
to  the  square  root  of  c,  which  is  approximately  pro- 
portional to  the  vertical  wavelength.  In  the  real 
atmosphere  the  meridional  heating  profile  q(y)  in  (12) 
does  not  adjust  to  different  vertical  scales  as  we  have 
assumed  in  our  calculations.  A  proper  treatment  of 
this  problem  is,  perhaps,  to  multiply  the  previously 
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calculated  response  by  the  projection 

j     q{y)expl-l3f/i2c)yy 
Jo 
P= 


exp[-/3>'V(2c)]2Jy 


(16) 


Here  we  could  specify  q{y)  as  we  did  for  the  vertical 
heating  profile,  although  there  appears  to  be  more 
observational  uncertainty.  For  the  purpose  of  the 
present  qualitative  discussion,  however,  this  is  not 
strictly  necessary.  We  know  that  the  latent  heating 
is  largely  confined  to  the  tropics,  so  that  the  meridional 
extent  of  q{y)  must  be  quite  limited.  Therefore,  the 
efTect  of  (16)  should  be  such  that  waves  whose  vertical 
scales  are  larger  than  that  corresponding  to  the  me- 
ridional scale  of  q{y)  will  suffer  from  having  a  large 
portion  of  their  north-south  domain  unsupported  by 
heating.  If  we  assume  that  the  meridional  scale  of  q{y) 
can  be  represented  by  an  e-folding  width  yd—  1500  km 
of  a  Gaussian  distribution  about  the  equator,  and 
that  6",  =  1.22X10-"  s-^  which  corresponds  to  r  =  3  K 
km-',  then  those  waves  with  a  tropospheric  vertical 
wavelength 

Zj>7r/35rH'd-~  14.2  km, 

or  a  stratospheric  wavelength 

L,>^.2  km, 

will  have  their  amplitude  reduced  by  (16).  The  reduc- 
tion will  be  more  severe  for  longer  vertical  wave- 
lengths. This  obviously  would  also  shorten  the  pre- 
viously calculated  most  favored  vertical  scales  of  the 
Kelvin  waves. 
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In  a  recent  paper  by  Kuo  (1975),  a  condition  was  waves  and  the  Kelvin  waves  in  the  tropics.  This  condi- 
used  in  a  CISK  parameterization  to  explain  the  time  tion,  as  stated  by  Kuo  on  p.  2231,  is  that  ''the  depletion 
and  length  scales  of  the  observed  mixed  Rossby-gravity      of  moisture  through  precipitation  must  be  replenished 
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by  evaporation  from  the  surface  during  a  complete 
period  of  the  disturbance."  Eq.  (15b)  in  his  paper 
specifies  this  condition  as  follows: 

r  (£+rate  of  moisture  convergence)  ^  maximum  total 
precipitation  given  by  deep  cumulus  convection 
during  a  period, 

where  r  is  the  period  of  the  wave  disturbances  and  E 
the  evaporation  rate.  On  p.  2234  he  derived  the  limiting 
period  tq  for  which  the  equal  sign  in  (15b)  holds.  This 
is  given  by  his  Eq.  (23) : 

Pt 

E 

where  Pt  is  the  average  maximum  total  precipitation 
over  a  period.  This  maximum  total  precipitation, 
according  to  the  argument  on  p.  2232,  is  always  pro- 
duced by  the  disturbance  "once  the  convective  insta- 
bility is  activated  by  the  large-scale  flow,"  because  the 
"convection  will  proceed  as  an  overturning  process  and 
produce  the  maximum  amount  of  precipitation  possible 
for  the  given  air  mass."  The  value  of  Pt  was  estimated 
by  subtracting  the  vapor  contents  of  the  outflow  column 
from  the  inflow  column  as  was  schematically  illustrated 
in  his  Fig.  3.  The  computed  ro  was  then  used  as  a  cutoff 
period  in  the  CISK  parameterization  [Eq.  (24)]  in  such 
a  way  that  the  heating  function  is  proportional  to  the 
factor  r/ro  if  r<ro.  The  implication  of  his  analysis  is 
that,  "disturbances  with  periods  <  xo  cannot  be  sup- 
ported fully  by  deep  cumulus  convection  initiated  by 
such  disturbances  because  the  moisture  content  of  the 
atmosphere  will  be  depleted  by  the  convective  process." 
The  above  arguments  are,  however,  inconsistent  with 
an  essential  part  of  the  CISK  parameterization  used  in 
the  paper,  namely  that  the  large-scale  heating  rate  is  pro- 
portional to  the  vertical  velocity  Wb  at  a  low  reference 
level  Zb  when  IFj>0.  In  this  parameterization  when  the 
wave  amplitude  at  an  early  stage  of  the  development  is 
small,  Wb  is  small  and  so  is  the  heating  rate,  which 
implies  that  the  precipitation  rate  is  small.  As  the 
waves  grow,  Wb  increases  and  the  precipitation  rate  in 
turn    increases.    Since    the    total   precipitation    of    a 


continuous  wave  chain  cannot  exceed  the  total  evapora 
tion  in  the  tropical  belt,  the  restraining  effect  is  sue 
that  the  precipitation  rate  must  not  exceed  the  evapora 
tion  rate,  or  3.75  g  cm"^  day~*  according  to  the  calcula 
tions  quoted  by  Kuo.  Therefore  the  total  moistur 
availability  in  the  tropics  should  constrain  the  max; 
mum  amplitude  the  waves  can  grow  to,  rather  than  th 
period  of  the  waves.  The  intensity  of  the  cumulus-seal 
convection  is  determined  by  the  large-scale  moistur 
convergence  below  Zb  in  the  CISK  theory  because  thi 
convergence  provides  the  high  6e  for  the  ascending  air 
For  more  sophisticated  parameterization  theories,  sucl 
as  that  of  Arakawa  and  Schubert  (1974)  which  is  i 
form  of  the  convective  adjustment,  the  same  incon 
sistency  would  occur.  In  these  theories  the  effect  o 
cumulus  convection  is  also  entirely  a  response  to  tht 
large-scale  forcing  and  therefore  its  intensity,  as  re 
fleeted  by  the  net  large-scale  heating,  should  be  con 
trolled  by  the  large-scale  flow. 

On  the  other  hand,  another  condition  mentioned  bj 
Kuo  [Eq.  (15a)3  is  a  valid  one,  and  may  be  quit( 
important  in  selecting  the  wave  scales.  This  conditior 
requires  "a  disturbance  of  sufficiently  long  duration  and 
high  intensity  able  to  lift  the  lower  atmosphere  above 
the  condensation  level,  thereby  activating  the  deep 
cumulus  convection."  In  other  words,  waves  with  short 
periods  are  unlikely  to  grow.  Kuo's  Eq.  (24),  which 
contains  a  cutoff  period  may,  in  a  way,  be  viewed  as  a 
crude  representation  of  this  condition  although  the 
precise  requirement  can  only  be  treated  as  a  nonlinear 
problem.  Based  on  this  condition  alone,  waves  whose 
Doppler-shifted  phase  speed  equals  zero  for  z^Zb  are 
most  favored  for  development.  This  seems  to  be  the 
case  in  the  tropical  western  Pacific  where  the  usuallyjj 
observed  warm-core  synoptic-scale  waves  have  a! 
typical  easterly  phase  speed  of  5-8  m  s~^,  quite  close  to 
the  low-level  mean  zonal  wind. 
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ABSTRACT 

In  this  paper  we  deal  with  the  interpretation  of  observed  oscillations  in  the  tropical  troposphere  and 
stratosphere  within  the  framework  of  the  equatorial  wave  theory.  A  dilHculty  with  this  problem  arises 
when  one  compares  the  short  vertical  wavelength  (or  equivalent  depth)  predicted  by  the  classical  theory 
and  the  observed  large  vertical  scales  associated  with  the  low  Doppler-shifted  frequencies  of  the  tropospheric 
oscillations.  In  this  analysis  it  is  shown  that  the  inclusion  of  simple  linear  damping,  justified  by  budget  studies 
which  revealed  the  important  role  of  cumulus  momentum  transport,  has  a  strong  influence  at  low  fre- 
quencies on  the  forced  equatorial  waves  and  results  in  two  types  of  dispersive  relationships.  The  first  type 
is  characteristic  of  the  regular  internal  gravity  waves  which  have  fast  phase  speeds  and  weak  vertical 
attenuation.  The  second  type  is  dominated  by  the  viscous  damping  time  scale  and  has  slow  phase  speeds 
and  strong  vertical  trapping.  The  theory  predicts  that  the  stratospheric  oscillations  may  be  identified  with 
the  first  type  and  the  tropospheric  oscillations  with  the  second.  In  the  case  of  Kelvin  waves  the  results 
can  be  used  to  explain  consistently  both  the  observed  stratospheric  Kelvin  waves  and  the  planetary-scale 
Kelvin-like  oscillations  in  the  troposphere  including  the  40-50  day  oscillation  and  the  monsoon  and  Walker 
circularions.  Possible  implications  with  respect  to  other  waves  in  the  tropics  are  also  discussed. 


1.  Introduction 

Based  on  spectral  and  cross-spectral  analyses  of 
long  time  series  (5-10  years)  of  tropical  station  data, 
Madden  and  Julian  (1971,  1972)  detected  a  global- 
scale  oscillation  in  the  tropical  troposphere  with  a 
periodicity  of  40-50  days.  A  degree  of  stationarity 
of  this  oscillation  in  time  was  also  established  by 
them  by  examining  the  limited  amount  of  station 
pressure  data  available  from  the  1890's.  Their  results 
indicate  that  this  oscillation  is  most  prominent  in  the 
zonal  velocity  component  and  the  surface  pressure, 
with  the  maximum  amplitude  situated  along  the 
equator.  There  are  some  time  and  spatial  variations 
within  a  40-50  day  cycle  as  depicted  in  Fig.  1  which 
is  reproduced  from  a  schematic  diagram  constructed 
by  Madden  and  Julian  (1972),  but  basically  the  oscil- 
lation propagates  eastward  at  a  fairly  steady  phase 
speed  and,  at  most  times,  exhibits  a  zonal  wave- 
number  1  structure.  Madden  and  Julian  (1972)  argued 
that  enhanced  deep  cumulus  convection  is  usually 
associated  with  large-scale  upward  motion  and  the 
oscillation  appears  as  two  vertical  circulation  cells 
along  the  equatorial  zonal  plane. 

The  above  structure  resembles  that  of  a  theoretical 
atmospheric  Kelvin  wave  (Hoi ton  and  Lindzen,  1968) 
in  several  ways:  the  direction  of  propagation,  the 
meridional  structure  and  the  absence  of  large  amplitude 
in   the  meridional  wind  component.   However,  for  a 


given  phase  speed  as  slow  as  that  observed,  Kelvin 
wave  theory  would  predict  a  fairly  short  vertical 
wavelength  which  does  not  appear  to  agree  with  that 
of  the  observed  oscillation.  This  discrepancy,  which 
has  been  a  major  difficulty  for  the  theoretical  under- 
standing of  the  40-50  day  oscillations,  has  also  been 
encountered  by  several  meteorologists  in  their  efforts 
to  interpret  various  t\-pes  of  large-scale  tropical  wave 
motions  using  the  equatorial  wave  theory.  In  a  linear 
analysis  Lindzen  (1967)  has  shown  that  for  synoptic 
(and  larger)  time  scales,  waves  on  an  equatorial 
beta-plane  generally  propagate  in  the  vertical  as  in- 
ternal gravity  waves  with  rather  short  vertical  wave- 
lengths. This  result  was  evident  in  a  study  of  the 
response  of  the  tropical  atmosphere  to  stationary 
heating  by  Webster  (1972,  1973).  In  a  two-level 
(750  and  250  mb)  numerical  model  he  found  that  the 
equatorial  response  in  an  easterly  basic  current  has 
a  horizontal  structure  that  is  characteristic  of  Kelvin 
waves.  However,  his  analytical  solution  for  the  case 
of  an  idealized  basic  flow  possesses  very  short  vertical 
wavelength  [~0(1  km)]  and  appears  to  suggest  that 
higher  vertical  resolution  is  required  for  the  nu- 
merical model. 

The  short  vertical  wavelength  predicted  by  the 
equatorial  wave  theory  arises  from  the  internal  gravity 
wave  character  of  these  waves.  Since  for  pure  gravity 
waves  the  phase  speed  c  is  given  by  c={gh)^,  where 
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Fig.  1.  Schematic  depiction  of  the  time  and  space  (zonal 
plane)  variations  of  the  circulation  cells  associated  with  the 
40-50  day  oscillation.  The  mean  pressure  disturbance  is  plotted 
at  the  bottom  of  each  chart  with  negative  anomalies  shaded. 
Regions  of  enhanced  large-scale  convection  are  indicated  sche- 
matically by  the  cumulus  and  cumulonimbus  clouds.  The  relative 
tropopause  height  is  indicated  at  the  top  of  each  chart.  (From 
Madden  and  Julian,  1972.) 


g  is  the  gravitational  constant  and  k  a  depth  scale 
for  the  motion  (usually  called  the  equivalent  depth), 
the  vertical  wavelength  tends  to  be  proportional  to 


the  phase  speed.  Within  the  range  of  large-scale 
tropical  motions  in  the  tropical  troposphere  the  phase 
speed  is  usually  slow  compared  to  the  external  gravity 
wave  phase  speed,  so  the  vertical  wavelength  is  con 
sequentially  small  compared  to  that  given  by  the 
scale  height.  The  observed  Kelvin  waves  (Wallace 
and  Kousky,  1968)  and  mixed  Rossby-gravity  wave 
(Yanai  and  Maruyama,  1966)  in  the  lower  strato- 
sphere have  rather  large  Doppler-shifted  phase  speeds 
(relative  to  the  mean  zonal  wind)  so  that  their  ob- 
served vertical  wavelengths  are  in  better  agreement 
with  the  theory.  In  Table  1  we  summarize  the  com- 
parison between  the  observed  vertical  wavelength^ 
and  that  given  by  the  equatorial  wave  theory  for 
three  major  zonal  wind  oscillations  in  the  tropics 
the  stratospheric  Kelvin  wave,  the  40-50  day  oscilla- 
tion and  the  planetary-scale  stationary  motion  which 
could  include  the  east-west  overturning  Walker  cir- 
culation along  the  equatorial  plane  (Bjerknes,  1969; 
Webster,  1973)  or  the  winter  monsoon  (Krishnamurti 
et  al.,  1973;  Webster,  1973).  The  theoretical  vertical 
wavelength  is  calculated  using  the  Doppler-shifted 
phase  speed  with  the  mean  zonal  flow  specified  by 
an  averaged  tropospheric  easterly  speed  of  —  5  m  s~"\ 
This  specification  is  based  on  the  evidence  that  all 
these  oscillations,  including  the  stratospheric  Kelvin 
wave,  are  generated  by  energy  sources  in  the 
troposphere. 

The  equatorial  wave  theory  used  in  Table  1  is 
basically  in  viscid  because  it  does  not  include  a  sig- 
nificant frictional  effect  in  the  free  atmosphere.  In  a 
numerical  diagnostic  study  of  the  seasonal-mean  cir- 
culation at  200  mb  during  the  northern  summer, 
Helton  and  Colton  (1972)  found  that  a  very  strong 
damping  term  (corresponding  to  linear  damping  time 
scale  of  ~  1  day)  is  needed  in  the  linearized  barotropic  j 
vorticity  equation  to  balance  the  generation  of  plane- 
tary-scale vorticity  by  horizontal  divergence.  Since 
the  planetary-scale  divergence  field  is  largely  a  result 
of  the  strong  condensation  heating  in  the  middle 
troposphere  due  to  the  vigorous  cumulus  convection 
associated  with  the  summer  monsoon,  Holton  and 
Colton  hypothesized  that  the  vertical  mixing  of 
vorticity  and  momentum  by  cumulus  convection 
could  account  for  this  strong  damping.  This  hypothesis 
is  also  consistent  with  several  budget  studies  of  syn- 
optic-scale, cumulus-related  tropical  motions  (Reed 
and  Recker,  1971;  Williams  and  Gray,  1973;  Reed 
and  Johnson,  1974;  Chu,  1976;  and  others),  and  the 
numerical  diagnostic  study  by  Chang  et  al.  (1975). 
Recently  Stevens  et  al.  (1977)  also  show  that  this 
strong  damping  mechanism  is  responsible  for  reducing 
the  theoretically  calculated  amplitude  of  t.ie  tropical 
waves  forced  by  heating  to  the  observed  magnitude. 
Thus  the  inclusion  of  this  mechanism  in  the  equatorial 
wave   theory  seems   to   be   necessary  for   the  proper 
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Table  1.  Comparison  of  observed  tropical  zonal  oscillations  and  the  inviscid  equatorial  wave  theory.  The  theoretical  vertical  wavelengths 
are  based  on  the  observed  Doppler-shifted  phase  speeds  which  assume  a  basic  zonal  wind  of  m=  —  5  m  s~'.  See  text  for  details. 
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interpretation  of  the  observed  tropical  motions  forced 
by  cumulus  heating. 

The  purpose  of  this  paper  is  to  study  the  effect 
of  the  strong  damping  mechanism  called  for  by  pre- 
vious budget  studies  on  tropical  wave  dynamics.  The 
simplest  form  of  linear  drag  will  be  used  to  represent 
this  effect  as  was  done  by  Holton  and  Colton  (1972). 
We  will  show  that  this  effect  in  some  cases  is  more 
than  just  an  expected  reduction  in  the  amplitude  of 
the  thermally  forced  waves  or  a  near-trivial  correction 
in  the  frequency  equation.  Moreover,  in  these  cases 
it  can  ;iccount  for  most  of  the  discrepancies  in  Table  1 
and  provide  a  consistent  theoretical  basis  for  the 
interpretation  of  the  tropical  oscillations  in  both  the 
stratosphere  and  the  troposphere. 

2.  Basic  equations 

Since  the  vertical  structure  equation  of  all  equatorial 
waves  is  identical  to  that  for  two-dimensional  gravity 
waves,  we  will  consider  only  the  Kelvin  wave  case 
for  the  purpose  of  simplicity.  Setting  the  meridional 
velocity  zero  we  may  write  the  linearized  zonal  mo- 
mentum, meridional  momentum,  hydrostatic,  thermo- 
dynamic energy  and  continuity  equations  on  an 
equatorial  beta-plane,  with  a  single  zonal  wave- 
number  k  and  constant  phase  speec  c  relative  to  the 
basic  zonal  flow  as 


—ikcu  =  —  ik<l> — Du, 

/3y«= , 

dy 


■ikcT-\-wT  = DT, 

Cp 


iku+e''^—{e-''"w)  =0, 

dz 


(4) 
(S) 


where  u,  w,  T,  <i>  and  Q  are  the  perturbation  zonal 
velocity,  vertical  velocity,  temiperature,  geopotential 
and  the  diabatic  heating  rate,  respectively;  if  is  a 
constant  scale  height,  F  the  static  stabiUty,  Cp  the 
specific  heat  at  constant  pressure,  R  the  gas  constant, 
^  the  meridional  gradient  of  the  vertical  component 
of  earth's  vorticity,  y  the  meridional  coordinate;  and 
z=—H\n{p/po)  is  the  vertical  coordinate  with  p  the 
pressure  and  po  a  reference  pressure.  The  damping 
effect  is  represented  by  the  linear  drag  coefi&cient  D  in 
the  zonal  momentum  equation.  A  Newtonian  cooling 
term  with  the  same  coefficient  is  included  in  the 
thermodynamic  energy  equation  in  order  to  facilitate 
the  analysis. 

Eqs.   (l)-(2)  give  the  meridional  structure  of  the 
Kelvin  waves : 


where 


M,<^a:exp(-/3y/2c), 


c=c+iD/k, 


(6) 


n)  and  the  condition  Re(c)  =  c>0  is  required  to  satisfy 
the  trapping  condition  on  an  equatorial  /3-plane.  If  a 
height-dependence  factor,  exp(z/2H),  is  separated  from 

(2)  the  perturbation  quantities,  Eqs.  (1)  and  (3)-(5)  may 
be  combined  to  form  a  single  equation  in  w : 


d<l>    RT 
dz      h' 


(3) 


dhu'  Q' 

hxv=-, 


(7) 


3t.- 
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The  parameter  X  is  a  measure  of  the  vertical  wave- 
number.  If  the  heating  function  Q'  is  given,  (7)  can 
be  solved  with  suitable  boundary  conditions.  In  our 
problem  the  following  boundary  conditions  are  used: 


w  = 


where 


0,  at     2  =  0  (9a) 

iCie'^'+Cie-'^',     at    z  =  Zt  (9b) 

Ci  =  rCz- 


Here  z<  is  the  height  of  tropopause  and  the  condi- 
tion (9b)  results  from  the  requirement  that  latent 
heating  vanishes  at  s^  The  parameter  r  is  a  reflection 
coefficient  which,  in  the  absence  of  vertical  wind 
shear,  is  given  by  the  model  specification  of  the  static 
stabihty  distribution.  The  solution  to  (7)  may  be 
written  in  the  form 


Az_i_  p — 1^2 


w{z)  = 


^giXzJ^g 


■w{z)  = 


A(l+r) 
sinAs 


smXs  — dz, 
gh 


z>. 


(10a) 


A(l+r) 


Q' 

gh 


re'^'+e 


—  1X2 


X 


smAz  — dz 
gh 


A(l+r) 

Z,>Z>Z<:,       (10b) 


w{z)=  — 


sinAz 

A(l+r)' 


Q' 


{re'^'+e-'^')—dz,  z.^z^O.   (10c) 
gh 


Here  Zt  is  the  height  of  cloud  top  and  Zc  the  height 
of  cloud  base.  Chang  (1976)  has  used  the  above 
equations  to  show  that,  in  the  case  of  D=0,  the 
vertical  energy  flux  of  Kelvin  waves  due  to  pressure 
work  is 

Si -2 


<j>'w' 


w 


(11) 


where  the  overbar  denotes  the  zonal  average  and 
(j)'  =  <j)exp(  —  z/2H).  Eqs.  (11)  and  (7)  indicate  that 
for  this  case  the  upward  forcing  of  Kelvin  waves  is 
most  efficient  for  long  waves.  Augmented  by  the  fact 
that  waves  with  a  vertical  scale  equal  to  twice  that 
of  forcing  are  most  efficiently  excited    (Green,   1965 


and  Lindzen,  1966),  Chang  (1976)  was  able  to  she 
that  the  Kelvin  wave  response  given  by  the  invisc 
theory  agrees  with  the  observed  stratospheric  Kelv 
waves. 

If  we  consider  downward  phase  propagation  (ui 
ward    energy   propagation)    only,    Re(A)>0.    For  . 
~7  km  and  vertical  wavelength  «88  km,   Eq.   (' 
may  be  approximated  by 


cX  =  SK 


(i: 


In  the  inviscid  case,  D=0  and  c=c  so  that  A  is  pu 
real.  In  the  case  of  D^^O,  c  =  Cr+iCi  where  Cr=c> 
and  Ci=D/k>0;  thus  A  becomes  complex,  i.e.,  X=, 
-fzAi.  Here  K  describes  the  exponential  increa, 
(A,>0)  or  decrease  (A,<0)  of  the  amplitude  wi| 
height.  Eq.  (12)  now  becomes  quadratic  and  has  tl^ 
following  roots : 

5^zt(5-4A,W)^ 

Crt   Tt= >0, 


'"•I. If 


^'l.ll'- 


2\r 


-<0. 


2ci 


The  subscripts  I  and  II  are  being  used  to  denote  tli' 
plus  and  minus  roots,  respectively.  It  is  obvious  tha 
Cri>Crii  and  |  A,i|  <  j  A,ii| ,  hence  the  first  solutio; 
(which  shall  be  called  Mode  I)  has  a  fast  horizonts 
phase  speed  and  can  propagate  more  efficiently  int 
higher  levels  while  the  second  solution  (Mode  II)  i 
slower  and  attenuates  faster  above  the  source  regior, 
In  addition,  Eqs.  (13)  and  (14)  indicate  that  the  wave! 
are  excitable  only  if 


S'w     kS'^ 


6* /     k^^\ 

A.<—    =—  , 

2c^\     2DI 


i 


so   that  the  effect  of  damping  causes  a  short- wavf 
cutoff  in   the  vertical  or  a  long-wave  cutoff  in 
horizontal. 

If  we  take  the  mviscid  limit  of  D-^,  Eqs.   (13) 
and  (14)  reduce  to 


Cri=5VAr 

A.i=0 

Crii=0 

A,„  =  -=o 


(16) 


(17)i 


Comparing  (16)  and  (12)  it  is  seen  that  Mode  I  cor- 
responds to  the  regular  inviscid  solution.  Substituting 
(17)  into  (l)-(5)  leads  to  a  trivial  solution  of  vanish- 
ing wave  amplitude.  Thus  Mode  I  may  be  called  the 
"regular  mode"  and  Mode  II  the  "viscous  mode,' 
because  the  latter  is  non-trivial  only  when  Z)?^0. 
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Fig.  2.  Real  Doppler-shifted  phase  speed  c,  as  a  function  of  vertical  wavelength  Lt  and  damping 
time.  For  the  convenience  of  application  to  other  types  of  equatorial  waves,  the  abscissa  can  also  be 
designated  by  (gA)r*,  where  h  is  the  equivalent  depth. 


Fig.  2  shows  the  real  phase  speed  Cr  as  a  function 
of  vertical  wavelength  Lz{  =  2T,\r)  and  the  viscous 
damping  time  scale  l/D.  For  £>  =  0  (damping  time 
scale  X )  the  dispersion  relationship  of  Mode  I  is  a 
straight  line  of  Ct"^  Lz,  typical  of  the  internal  gravity- 
waves.  In  addition,  the  vertical  coordinate  (cr=0) 
represents  the  trivial  solution  of  Mode  II.  For  Z)r^O 
the  two  modes  are  represented  by  the  quadratic 
curves.  The  effect  of  damping  on  Mode  I  is  a  rela- 
tively small  modification  with  Lz  increasing  as  Cr  in- 
creases. The  modification  for  Mode  II  is  more  sub- 
stantial and  Lz  generally  decreases  as  Ct  increases. 
The  two  modes  are  separated  by  the  line  Cr=c,.  For 
Mode  I  CT>Ci,  i.e.,  the  inertial  time  scale  is  faster 
than  the  damping  time  scale.  The  reverse  is  true  for 
Mode  II. 

The  imaginary  vertical  wavenumber  X,  is  shown  in 
Fig.  3  as  a  function  of  vertical  wavelength  and  damping 
time  scale.  Here  the  magnitude  of  X,  decreases  with 
increasing  Lz  for  Mode  I  and  increases  with  Lz  for 
Mode  II.  The  Cr=Ct  curve  again  provides  the  parti- 
tion of  the  two  modes. 

From  Figs.  2  and  3  it  can  be  seen  that,  if  a  certain 
vertical  wavelength  is  most  efficiently  excited  by  a 
given  forcing  function,  both  modes  will  be  excited 
for  Dj^O.  The  regular  mode  will  have  fast  zonal  phase 
speed  and  less  vertical  attenuation,  while  the  viscous 
mode  will  have  slow  zonal  phase  speed  and  is  more 
trapped  in  the  vertical  away  from  the  source  region. 


3.  Response  of  specified  heating 

The  Kelvin  wave  response  of  the  tropical  atmo- 
sphere can  be  calculated  if  the  vertical  heating  func- 
tion Q'  in  (7)  is  known.  As  in  Chang  (1976)  we  will 
use  the  vertical  heating  profile 


(3'cc5^(y)  sinirl 


<2'=o, 


(?)■  ^ 


(^S^Cc 


2>St  or  s<z<; 


where  the  function  gCy),  which  specified  the  meridio- 
nal distribution  of  heating,  is  assumed  to  take  the 
form  of  (6).  Solution  (10)  will  be  evaluated  by  speci- 
fying the  vertical  scale  of  the  heating  Az=St— 2c 
=  14  km,  Z)  =  2.3XlO-«  s-'i  (damping  time  =  5  days) 
and  5,=35t,  where  5,  and  5t  are  the  static  stability 
for  the  stratosphere  and  troposphere,  respectively. 

Since  Mode  I  is  less  trapped  in  the  vertical  we  ex- 
pect that  it  may  be  more  prominent  away  from  the 
source  region.  Using  (10a)  we  display  the  vertical 
wave  energy  flux  due  to  pressure  work,  ^'w' ,  at  Zt  as 
a  function  of  vertical  wavelength  Lz  and  zonal  wave- 
number  s{  =  ka,  where  a  is  the  radius  of  the  earth), 
in  Fig.  4.  The  cross-hatched  area  in  the  lower  part 
of  the  diagram  indicates  the  short  vertical  wavelength 
cutoff  required  by  (15).  The  maximum  response  occurs 
at  zonal  wavenumber  1-2  and  vertical  wavelength 
30-35  km.  As  expected,  this  result  is  similar  to  that 
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Fig.  3.  Imaginary  vertical  wave-number  X,  as  a  function  of  vertical  wavelength  and  damping  time. 


obtained  by  Chang  (1976)  except  that  the  vertical 
wavelength  is  slightly  longer  and  wavenumber  2  be- 
comes comparable  in  importance  to  wavenumber  1. 
This  vertical  wavelength  range,  according  to  Fig.  2, 
implies  a  Doppler-shifted  phase  speed  ~50  m  s~',  in 
reasonable  agreement  with  the  observed  stratospheric 
Kelvin  waves. 

For  the  highly  trapped  Mode  II  the  vertically- 
integrated  tropospheric  kinetic  energy  is  considered 
an  appropriate  measurement  for  the  response  spectrum. 


Fig.  5  shows  this  quantity  as  a  function  of  L,  and  s. 
Here  the  zonal  scale  selection  clearly  favors  the  lowest 
zonal  wavenumber,  but  all  vertical  wavelengths 
beyond  the  short-wave  cutoff  have  about  the  same 
response.  Thus  the  most  prominent  viscous  mode  ex- 
cited should  be  of  zonal  wavenumber  1  and  vertical 
wavelength  >  17  km.  For  vertical  wavelength  "-17  km 
the  Doppler-shifted  phase  speed  is  ~10  m  s~'  accord- 
ing to  Fig.  2.  These  values  agree  quite  well  with  the 
40-50  day  oscillation  observed  by  Madden  and  Julian 


*■>*'  for  MODE  I 


Fig.  4.  Vertical  wave  energy  flux  due  to  pressure  work  at  the  tropopause  level, 
<t>'w' ,  as  a  function  of  vertical  wavelength  Li  and  zonal  wavenumber  s  for  Mode  I. 
Units  are  arbitrary. 
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K.E.  Response  for  MODE  II 


Fig.  5.  Tropospheric  kinetic  energy  response  for  Mode  II  as  a  function  of  vertical 
wavelength  and  zonal  wavenumber.  Units  are  arbitrary. 


(1972).  However,  if  the  forcing  function  has  a  white 
noise  distribution  in  frequency  the  maximum  Mode  II 
response  is  still  somewhat  smaller  than  the  maximum 
Mode  I  response.  The  absence  of  detectable  10-15 
day  Kelvin  waves  in  the  troposphere  therefore  indi- 
cates that  the  thermal  forcing  is  of  a  red  noise  nature. 
However,  this  requirement  of  red  noise  distribution  is 
not  so  much  for  the  explanation  of  the  stratospheric 
Kelvin  waves,  as  suggested  by  Holton  (1973);  rather, 
it  is  for  the  explanation  of  the  dominance  of  low- 
frequency  oscillations  in  the  troposphere.  The  10-15 
day  Kelvin  waves  in  this  case  are  of  relatively  higher 
frequencies  and  can  become  prominent  only  above  the 
troposphere  where  the  lower  frequency  viscous  mode 
is  greatly  attenuated. 

The  above  results  are  summarized  in  Table  2  within 
the  "non-resonant  selection"  categorv'.  The  stationary 
waves  related  to  the  Walker  circulation  or  the  winter 
monsoon  may,  on  the  other  hand,  be  viewed  as  due 
to  resonant  excitation  by  the  stationary  heat  sources, 
because  the  forcing  function  has  a  pre-selected  zero 
ground  phase  speed.  By  assuming  a  zonal  mean  wind 


of  —  5  m  s~^  in  the  source  region  the  Doppler-shifted 
phase  speed  of  this  resonance  is  5  m  s~^  Referring 
back  to  Fig.  2,  response  is  possible  only  in  the  Mode  II 
regime  with  a  vertical  wavelength  of  ~27  km  for  the 
case  of  5-day  damping  time.  This  vertical  wavelength 
agrees  much  better  with  the  observed  monsoon  or 
Walker  circulation  than  Webster's  (1972)  inviscid 
calculation  of  ~1  km.  In  addition.  Fig.  5  still  indicates 
that  zonal  wavenumber  1  should  be  clearly  preferred 
over  shorter  waves  if  the  heat  source  has  a  white 
noise  distribution  in  the  zonal  scales.  In  the  real 
atmosphere  the  differential  heating  due  to  land-sea 
differences,  as  in  the  case  of  monsoon,  or  due  to  sea- 
surface  temperature  variations,  as  in  the  case  of 
Walker  circulation,  is  usually  confined  to  lower  zonal 
wavenumbers ;  so  it  is  not  surprising  to  observe  the 
large  zonal  scale  of  these  circulations.  On  the  other 
hand,  Egger  (1976)  reported  that  the  response  to  the 
imposed  heating  anomalies  in  a  linear  model  almost 
always  favors  the  lowest  wavenumbers  regardless  of 
the  zonal  scale  of  the  anomaly  imposed  on  the  model. 
His  experimental  result  is  entirely  understandable  in 


Table  2.  The  most  efficiently  e.xcited  Kelvin  waves  predicted  by  the  viscous  theory. 


Vertical 

wavelength 

(km) 

Zonal 
wavenumber 

Doppler-shifted 

phase  speed 

(m  s-') 

Ground 

phase  speed 

(m  s-') 

Resonant  selection 
Non-resonant  selection 

Mode    I 
Mode  II 
Mode    I 
Mode  II 

31-35 

(not  excitable) 

27 

1-2 
1 

1 

50-58 
15 

5 

45-53 
0 

Damping  time  =  5  days. 
Vertical  scale  of  heating  =  14  km. 
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view  of  our  result.  This  low-wavenumber  preference 
for  the  thermally  forced  Kelvin  wave  response  is 
related  to  (11)  and  the  explanation  has  been  given 
by  Webster  (1973)  and  Chang  (1976). 

4.  Further     comparison     with     the     40-50     day 
oscillation 

As  mentioned  previously,  the  two  solutions  of  (12) 
are  separated  by  the  curve  c,=c,.  Since  for  Mode  I 
the  inertial  time  scale  is  faster  than  the  damping  time 
scale,  it  follows  that  the  main  balance  in  Eq.  (1)  is 
between  the  first  two  terms,  i.e., 

—ikcu^ik<f>. 

In  this  case  u  and  </>  are  in-phase  and  the  wave  is 
characteristic  of  the  internal  gravity  waves.  For 
Mode  II  the  damping  time  scale  is  faster  so  that  the 
main  balance  is  between  the  last  two  terms,  i.e., 

0=  —  ik<t>—Du. 

Thus  u  and  <f>  have  a  quarter-cycle  phase  difference 
which  is  characteristic  of  a  highly  viscous  wave  motion. 
If  we  examine  the  schematic  diagram  of  the  40-50  day 
oscillation  (Fig.  1)  depicted  by  Madden  and  Julian, 
It  can  be  seen  that,  for  most  of  the  categories,  u  and  <t> 
have  a  phase  relationship  between  in-phase  and 
quarter-cycle  difference.  The  important  effect  of 
damping  is  thus  clearly  indicated. 

If  the  vertical  scale  is  large,  the  main  balance  in 
the  thermodv-namic  energy  equation  is  between  the 
adiabatic  cooling  and  the  diabatic  heating.  The  residue 
is  approximately  balanced  by  the  advective  term  for 
Mode  I,  i.e., 

—  ikcT^ 'd.'V 

and  T,  w  will  have  a  quarter-cycle  phase  difference. 
For  Mode  II  this  residue  is  mainly  balanced  by  the 


damping  term,  i.e.. 


wV '=iDT. 


In  this  case  T,  w  will  be  in-phase  if  diabatic  heating 
exceeds  adiabatic  cooling  and  half-cycle  out-of-phase 
if  adiabatic  coohng  dominates  over  diabatic  heating 
In  Madden  and  Julian's   (1972)  observational  study 
they   found    that    the   surface   pressure   is   half-cycle 
out-of-phase  with  the  300  mb  temperature  and  in-phase 
with  the  100  mb  temperature.  Since  vertical  motion 
IS  close  to  in-phase  with  surface  pressure  as  shown 
by  Fig.  1,  It  may  be  inferred  that  T  and  w  are  in-phase 
at  300  mb  and  half-cycle  out-of-phase  at  100  mb.  The 
former  level  is  in  the  upper  troposphere  where  diabatic 
heating  must  dominate   for   the   motion   to   be   self- 
sustaining,  while  the  latter  level  is  near  the  tropopause 
where  the  diabatic  heating  vanishes  but  the  vertical 
motion  may  overshoot  through   the  cloud   top.  The 
thermal  structure  of  the  40-50  day  oscillation   thus 
also  appears  to  be  consistent  with  the  viscous  theory. 
When  damping  is  included,  the  meridional  structure 
of  the  Kelvin  waves  as  expressed  by  (6)  contains  a  tilt 
of  the  wave  axis  because  c  is  complex,  i.e., 

-=*(-- ^j-p(,---yj.    (18)  ^ 

Here  the  first  part  of  the  right-hand  side  gives  a 
Gaussian  distribution  in  the  meridional  direction  with 
a  Gaussian  folding  with  yc  given  by 

'2S' 


'°=fci;)''&''"-)" 


T^ 


'""  90*W  5; ' 

^ro::.^'enlZ^::::^n^^^^  f-  apprcx^atelv  the  36-50  day  penod  range  of 
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FiG.  7.  Comparison  of  the  meridional  phase  tilt  of  the  40-50  day  oscillation  at 
140°E,  170°N,  observed  by  Madden  and  Julian  (1972),  and  the  value  of  the  most 
efi&cient  Mode  II  response  predicted  by  the  viscous  theory. 


This  meridional  trapping  scale  is  identical  for  both 
the  regular  and  the  viscous  modes  for  a  given  vertical 
wavelength.  In  the  case  of  the  most  efficient  Mode  II 
response  discussed  in  the  previous  section,  \yG\ 
=  1500  km  and  agrees  quite  well  with  the  observed 
40-50  day  oscillation.  The  second  part  of  the  right- 
hand  side  of  (18)  gives  rise  to  a  periodic  variation  of 
the  waves  in  y,  and  therefore  a  meridional  tilt  of  the 
wave  axis.  It  also  implies  a  poleward  phase  propaga- 
tion for  the   Kelvin  waves  because  6r>0,  c,>0  and 

\r>0. 

A  strong  meridional  tilt  was  found  by  Madden  and 
Julian  for  the  40-50  day  oscillation  using  cross-spectral 
analysis  of  the  surface  pressure  data.  Their  phase 
diagram,  constructed  by  using  the  Canton  Island  series 
as  the  base  series,  is  reproduced  here  in  Fig.  6.  Al- 
though the  meridional  phase  tilt  seems  to  vary  in 
different  geographical  locations,  near  Canton  Island 
and  near  the  equator  the  tilt  is  primarily  NW-SE 
north  of  the  equator  and  NE-SW  south  of  the  equator. 
For  eastward  propagating  waves  this  indicates  a  pole- 
ward phase  propagation.  In  Fig.  7  the  phase  tilt  of 
the  40-50  day  oscillation  deduced  from  Fig.  6  at  two 
different  longitudes,  170°W  (near  Canton  Island)  and 
140°E  (active  convection  area  in  Indonesia),  are 
plotted  along  with  that  given  by  the  viscous  theory. 
Within  the  limitation  of  our  simple  theor\'  and  the 
approximations  invoked,  the  theoretical  result  can  be 
considered  as  fairly  good  in  simulating  the  observed 
phase  tilts. 


5.  Concluding  remarks 

We  have  shown  that  the  inclusion  of  a  damping 
term  can  result  in  two  types  of  dispersive  relationships 
for  the  forced  Kelvin  waves.  The  first  t>'pe  (Mode  I) 
is  characteristic  of  the  regular  internal  gravity  waves 
which  propagate  rapidly  in  the  zonal  direction  and 
attenuate  slowly  in  the  vertical.  The  second  type 
(Mode  II)  is  dominated  by  the  viscous  damping  time 
scale  and  has  slower  zonal  phase  speed  and  faster 
vertical  attenuation.  This  viscous  mode  reduces  to  a 
trivial  solution  in  the  inviscid  limit  so  its  existence 
is  supported  only  by  damping.  Assuming  that  the 
main  thermal  sources  in  the  tropics  are  confined  to 
the  troposphere  and  have  a  red  noise  distribution 
which  enhances  the  response  of  the  viscous  mode, 
this  theor\'  predicts  that  the  stratospheric  waves 
behave  like  the  regular  mode  and  the  tropospheric 
waves  behave  like  the  viscous  mode.  Our  theory  in- 
cludes features  that  are  consistent  with  several  results 
found  by  previous  theoretical  or  numerical  studies  of 
the  equatorial  waves  that  incorporate  the  damping 
mechanism.  Lindzen  (1971)  and  Holton  and  Lindzen 
(1972)  found  that  as  a  wave  is  propagating  in  weak 
vertical  wind  shear,  the  efficiency  of  damping  would 
change  inversely  according  to  the  change  of  the 
Doppler-shifted  phase  speed.  This  is  consistent  with 
our  result  that  the  slower  moving  waves  are  more 
trapped  in  the  vertical.  In  a  numerical  diagnostic 
model  of  the  Kelvin  waves,  Holton  (1973)  observed 
a  maximum  response  in  the  tropospheric  wavenum- 
ber  1  spectrum  with  period  >30  days.  This  response 
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seems  to  correspond  to  our  viscous  mode.  Very  recently 
Stevens  et  al.  (1977)  reported  that  with  a  parameter- 
ization of  the  cumulus  damping,  the  tropical  waves 
have  very  little  vertical  structure.  This  is  again  ex- 
plainable by  our  theory  because  the  vertical  wave- 
length is  increased  as  a  result  of  damping. 

In  real  atmosphere  the  damping  due  to  cumulus 
convection  and  other  processes  is  much  more  com- 
plicated than  the  simple  linear  drag  term  used  in  our 
analysis.  However,  the  agreement  with  observed  zonal 
wind  oscillations  such  as  the  stratospheric  Kelvin 
waves,  the  40-50  day  oscillations,  and  the  stationary 
monsoon  and  Walker  circulations  strongly  indicate 
that  our  theoretical  interpretation  is  quite  relevant 
to  these  oscillations.  In  addition,  although  the  present 
analysis  is  limited  to  Kelvin  waves  only,  the  vertical 
structure  equation  applies  to  all  equatorial  waves. 
Hence  it  would  be  quite  possible  that  the  general 
result  also  applies  to  other  thermally  controlled  tropical 
wave  disturbances.  This  is  particularly  interesting  in 
view  of  the  large  vertical  scale  and  slow  Doppler- 
shifted  phase  speed  usually  observed  in  the  tropical 
tropospheric  waves  such  as  the  easterly  waves.  In  a 
discussion  of  the  present  results,  Wallace  (personal 
communication)  suggested  that  they  would  also  ex- 
plain the  seemingly  non-dispersive  westward  propaga- 
tion of  the  synoptic-scale  cloud  clusters  observed  in 
the  time  sections  of  tropical  satellite  pictures  (Chang, 
1970).  In  these  sections  the  propagating  cloud  patterns 
have  quite  uniform  phase  speeds  but  the  apparent 
wavelength  varies  over  a  large  range.  This  non-dis- 
persive character  led  Wallace  (1972)  to  propose  that, 
despite  the  strong  vertical  and  longitudinal  variations 
of  the  basic  flow,  they  are  more  likely  a  result  of 
advection  by  basic  flow  than  a  manifestation  of  some 
kind  of  Rossby  waves.  He  argued  that  the  difficulty 
of  ascertaining  a  steering  level  is  not  as  severe  as  that 
of  explaining  why  the  beta-effect  does  not  produce 
dispersive  phase  speeds.  However,  to  interpret  these 
patterns  as  being  passive  cloud  masses  advected  by 
the  basic  flow  is  also  unsatisfactory  due  to  the  fact 
that  they  are  frequently  observed  to  be  associated 
with  active  convections  and  that  they  sometimes 
decay  and  regenerate  during  the  course  of  propaga- 
tion. These  difficulties  are  partially  alleviated  in  the 
framework  of  the  present  theory  because  the  viscous 
mode  should  have  a  much  smaller  Doppler-shifted 
phase  speed  and  would  appear  as  if  they  are  advected 
by  the  basic  easterly  current  at  some  levels. 
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Some  Theoretical  Problems  of  the 
Planetary-Scale  Monsoons 
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Siimimirx  -  Some  important  theoretical  problems  of  the  planetary-scale  monsoons  which  have  arisen 
from  recent  advances  of  observational  studies  are  reviewed.  These  include:  ( I )  the  requirement  ot"a  strong 
damping  mechanism  in  the  planetary  scale  vorticity  budget  of  summer  monsoon  and  a  similar  but  weaker 
requirement  tor  the  wimcr  monsoon:  (2)  the  localized  baroiropic  instability  of  the  summer  monsoon 
which  IS  a  result  of  the  strong  zonal  asymmetry  of  the  planetary-scale  flow  and  causes  significant  nonlinear 
energy  conversions;  and  (3)  the  oscillations  of  the  planetary-scale  monsoons.  It  is  pointed  out  that  these 
problems  are  inter-related  and  their  understanding  is  also  important  for  the  proper  simulation  of  other 
scales  of  motion  of  the  monsoon  circulation. 

Key  words:  Monsoon    Oscillation,  \onicity  budget:  Baroiropic  instability. 


/.    I/ur/xluclion 

The  monsoon  circulation  o\'  both  the  northern  summer  and  the  northern  winter 
(hereafter  the  word  "northern"  will  be  neglected)  are  characterized  by  large  planetary 
scale  quasi-stationar\  waves  in  the  subtropics  and  tropics.  As  have  been  shown  by 
Krishn.amlrti  (1971a),  Krishnami  rti.  Kana.mitsi  .  K(jss  and  Lee  (1973b)  and 
others,  these  waves  are  a  result  o\'  the  thermally  direct  overturnings  induced  by  the 
differential  heating  associated  with  the  land-sea  distribution.  The  major  component 
of  the  heating  is  the  convective  latent  heat  release  which  is  most  intense  over  the 
Tibetan  and  Mexican  highlands  during  summer  and  the  maritime  continent  of 
the  Malaysia-Indonesia  region  during  winter.  Mainly  due  to  the  availability  of  the 
commercial  aircraft  data,  the  principal  features  of  these  zonal  asymmetries  have  been 
best  depicted  by  the  wind  analyses  near  200  mb  by  Krishna.ml  RTiand  R()GERS(I970). 
Krishnamlrti  et  al.  (1973b)  and  Sadler  (1975).  During  summer  these  features 
include  a  strong  anticyclone  centered  over  the  Tibetan  plaleau.  a  weaker  anticyclone 
over  Mexico  and  elongated  troughs  oriented  northeast-southwest  in  the  m.id-Pacific 
and  mid-.Atlantic.  These  are  all  well  presented  in  Fig.  1  which  shows  the  time-mean 
streamfunclion  at  200  mb  for  June  .August  1967  as  determined  by  Krishnami  rti 
( 1971a).  During  winter  a  large  streamtunction  high  over  the  equatorial  South  China 
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Figure  1 
Summer  mean  (June -August  !967)  ZOO-mb  slreamfunction,  contour  interval  is  50   '    10*^  m-  sec"  '.  (From 

Krishnamurti.  1971a.) 


Figure  2 

Winter  mean  (  December -February  1971 )  20()-mb  streamtunction.  contour  interval  is  50    ■    10^  m'sec"'. 

(From  Krishnami-rti  £'/ tv/.  1973b.) 

Sea  and  western  Pacific  replaces  the  Tibetan  high  to  become  the  most  pronounced 
planetary  scale  feature,  with  two  weaker  highs  appearing  over  the  equatorial  eastern 
Africa  and  the  northeastern  South  America  and  three  oceanic  troughs  extending 
from  the  equator  to  the  southern  subtropics.  These  features  are  shown  in  Fig.  2  which 
is  the  200  mb  time-mean  streamfunction  for  the  winter  of  1971  determined  by 
Krishnamurti  et  al.  (1973b). 

A  number  of  interesting  theoretical  problems  have  arisen  from  the  observed 
monsoon  circulations.  In  this  paper  we  will  review  some  of  the  problems  that  are 
related  to  the  planetary-scale  monsoon.  They  include: 

( 1 )  The  requirement  oi  a  strong  damping  mechanism  in  the  planetary  scale 
vorticitv  budget. 
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(2)  The   localized    barotropic    instability   associated    with    the   planclary-scalc 
summer  monsoon,  and 

(3)  Oscillations  of  the  planetary-scale  monsoons. 


2.    The  requirement  of  a  strorui  planelary  seale  damftinq  fuedianisni 

In  a  linear,  one-level  numerical  model,  Holton  and  Coljon  (1472)  used  the 
barotropic  vorticity  equation  to  diagnose  the  budget  of  the  seasonal-mean  planetary 
scale  vorticity  at  200  mb  during  summer.  They  used  the  200  mb  mean  divergence 
observed  by  Krishnamurti  (1971a,  b)  as  the  steady  forcing  function  with  a  mean 
zonal  flow  also  specified  by  the  observed  values.  Figure  3  shows  the  velocity  potential 


Figure  ? 
Summer  mean  (June-August  1967)  200-mb  velocity  potential,  contour  interval  is  10  x   |0'  m"  sec^  '.  The 
tew  streamlmes  sketched  indicate  the  direction  ot  the  mean  divergence  motions.  (From  Krishnamurti. 

1971a.) 


computed  by  Krishnamurti  (1971a)  which  expresses  the  observed  mean  divergence 
field.  Comparing  this  figure  with  Fig.  1,  it  is  evident  that  the  time-mean  vorticity  and 
divergence  fields  are  approximately  half-cycle  out  of  phase.  However,  the  resultant 
steady-state  planetary  scale  vorticity  field  calculated  by  Holton  and  Colton  is 
about  one-quarter  cycle  out  of  phase  from  the  divergence  field,  unless  a  very  strong 
damping  term  is  added  to  the  mean  vorticity  equation.  This  difficulty  arises  because, 
in  the  absence  of  strong  damping,  the  main  balance  in  the  vorticity  equation  is  be- 
tween the  generation  by  horizontal  divergence  and  the  advection  by  mean  flow.  This 
near-balance  requires  that  the  vorticity  be  one-quarter  cycle  out  of  phase  from  the 
divergence.  The  very  efficient  generation  by  the  strong  divergence  also  causes  the 
amplitudes  of  the  long  waves  to  be  significantly  larger  than  those  observed.  Figure  4 
(from  Holton  and  Colton.  1972)  compares  the  diagnosed  and  observed  amplitudes 
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and  phases  of  the  streamfunction  for  zonal  wavenumber  oneasa  function  of  latitude. 
When  the  linear  damping  coefficient  is  increased  from  a  modestly  strong  value  of 
5  X  10~^sec~'  to  a  very  strong  value  of  1.5  x  10"-  sec"',  the  discrepancy  dis- 
appears and  the  phase  and  amplitude  of  the  long  waves  become  in  close  agreement 
with  observations. 
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Figure  4 

Amplitude  in  units  of  10*"  m^  sec"  '  (left-hand  curves)  and  phase  (right-hand  curves)  of  the  streamfunction 

for  zonal  wavenumber  1  as  a  function  of  latitude.  Solid  curve  is  observed  streamfunction;  short  and  long 

dashes  are  the  solutions  computed  by  Holton  and  Colton  (1972)  for  damping  rate  of  1.5  *   10" '  sec" ' 

and  5   '   10"^  sec" '.  respectively.  (From  Holton  and  Colton,  1972.) 

HoLTO.N  and  Colton's  (1972)  findings  suggest  that  a  very  strong  sink  with  a 
damping  time  scale  faster  than  1  day  is  required  to  balance  the  linear  vorticity  budget 
of  the  seasonal-mean  200  mb  planetary-scale  tlow  dunng  summer.  This  damping  is 
much  stronger  than  the  normal  dissipation  process  in  the  upper  troposphere  whose 
damping  time  is  on  the  order  of  weeks.  In  looking  for  an  explanation  of  this  excessive 
damping  requirement,  they  suggested  that  two  possible  mechanisms  not  included  in 
their  model  may  be  responsible  for  the  strong  damping.  The  first  is  the  vertical  mixing 
of  vorticity  by  cumulonimbus  convection  which  is  very  strong  over  the  Asian  conti- 
nent and  Mexican  highland  during  the  northern  summer.  This  mechanism  has  been 
suggested  by  Riehl  and  Pearce  (1968).  Reed  and  Recker  (1971),  Wallace  (1971). 
WiLLiA.vis  and  Gray  (1973)  and  others  to  balance  the  vorticity  budget  of  the  synoptic- 
scale  easterly  waves  in  ihe  tropical  troposphere.  However,  Holton  and  Colton 
(1972)  felt  that  it  alone  is  still  insufficient  to  account  for  the  required  fast  damping. 
They  then  proposed  the  second  mechanism  which  is  related  to  the  day-to-day 
fluctuations  of  the  monsoon  vorticity  and  divergence  fields  observed  by  Krishnamur- 
Ti  ( 1971b).  Considering  deviations  from  the  time  mean,  the  complete  divergence  term 
in  the  time-mean  vorticity  equation  is  <(C  ^  f)^^  =  (<C>  +  /)<(^>  +  <dC'>  where  C 
and  /are  the  relative  and  planetary  vorticities,  respectively;  <)  is  divergence  and  the 
primes  are  deviation  from  the  lime  mean  <   ).  The  last  term,  which  was  excluded  in 
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HoLTON  and  Colton's  linear  calculations,  tends  to  be  negative  according  to  the 
observed  out-of-phase  correlation  of  the  fluctuating  vorticity  and  divergence  fields. 
This  term  therefore  may  counteract  the  vorticity  generation  by  the  time-mean 
divergence.  Holton  and  Colton  then  speculated  that  the  combination  of  this 
fluctuation  correlation  and  the  cumulus  mixing  mechanisms  could  be  responsible  for 
the  less  than  one  day  damping  time. 

The  requirement  of  a  strong  vorticity  sink  has  also  been  confirmed  by  other 
numerical  investigations  that  do  not  include  the  strong  damping.  As  a  follow  up  of 
the  work  of  Holton  and  Colton  (1972).  Colton  ( 1973)  used  a  semi-spectral  repre- 
sentation of  the  one-level  nonlinear  barotropic  vorticity  equation  to  study  the  scale 
interactions  of  the  200  mb  summer  monsoon  flow  forced  by  the  observed  mean 
divergence.  As  a  consequence  of  the  inclusion  oi"  only  modest  damping  (time  scale 
5  days),  the  resultant  planetary-scale  waves  have  amplitudes  considerably  larger  than 
the  seasonal  mean  values.  Westward  phase  shifts  of  50°  for  the  Tibetan  high  and  30° 
for  the  mid-Pacific  trough,  in  comparison  with  the  observed  positions,  are  also 
produced  by  the  model  (Fig.  5).  Abbott  ( 1973)  developed  a  3-level  nonlinear  quasi- 
geostrophic  hemispheric  model  to  simulate  the  time-mean  climatology  of  the  summer 


Figure  5 

Streamlanction  of  the  forced  motion  averaged  over  days  15  through  -W)  of  the  numerical  integration  by 

Colton  (197.1).  Contour  interval  IS  50    -    lO^m'sec    '. 


monsoon.  To  drive  the  monsoon  circulation,  he  specified  a  steady  heating  distribution 
which  was  found  by  inverting  the  governing  equations  with  the  vorticity  field  giving 
by  the  observed  time-mean  values.  This  heating  distribution  (Fig.  6)  in  turn  success- 
fully simulates  the  prmcipal  features  of  the  planetary-scale  flow  as  well  as  the  nonlinear 
energy  exchanges.  However,  it  is  about  a  quarter  cycle  out  of  phase  from  that  implied 
by  the  observed  mean  divergence  field  and  the  satellite  cloud  data.  This  apparently  is 
also  due  to  the  absence  of  a  strong  damping  mechanism  in  the  model  which  forces  a 
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Figure  6 
Time-independent  horizontal  distribution  of  the  diabatic  heating  rate  at  350  mb  as  computed  by  Abbot 

(1973).  Isopieths  are  at  5   <   10''' m- sec  " '. 

balance  between  the  advection  and  generation  of  vorticity.  In  a  3-component  semi- 
spectral  primitive  equation  model,  Murakami  (1974)  also  used  a  specified  heating 
function  to  study  the  planetary  scale  monsoon  flow  during  summer.  Without  the 
inclusion  of  the  strong  dampmg  and  allowing  the  specified  heating  to  fiuctuate  with 
time,  he  was  able  to  obtain  the  proper  amplitude  and  phase  of  the  long  waves. 
However,  the  resultant  divergence  field  in  his  model  is  somewhat  different  from  that 
observed  by  Krishnamlrti  (1971a).  The  greatest  discrepancy  occurs  at  200  mb. 
where  Krishnamlrti  ( 1971a)  found  very  strong  divergence  in  the  Tibetan  high  area 
during  the  1967  summer,  while  in  Murakami's  ( 1974)  model  this  level  is  almost  non- 
divergent.  The  model  also  produces  a  different  horizontal  distribution  of  long  wave 
divergence  compared  to  the  observed  1967  distribution,  with  the  Tibetan  divergence 
center  shifted  slightly  to  the  east  and  the  mid-Atlantic  convergence  center  shifted  to 
the  longitude  of  central  Africa.  Thus  the  magnitude  of  divergence  forcing  is  much 
reduced  and  a  phase  shift  between  divergence  and  anticyclonic  vorticity  also  results' ). 
In  addition,  the  advective  effect  is  not  properly  treated  because  the  zonal  mean  How 
is  not  allowed  to  interact  with  the  three  interacting  long  waves.  Hence  it  is  unclear 
what  is  the  main  mechanism  that  produces  the  correct  amplitude  and  phase  of  the 
vorticity  field. 

Looking  for  alternative  sources  for  the  strong  damping,  Barcilon  and  C(x:)per 
(1974)  hypothesized  that  a  thermal  boundary  layer  effect  could  operate  in  the 
atmosphere  during  summer  monsoon.  In  their  simple  analytical  model  a  sinusoidal 
temperature  variation  is  applied  at  the  lower  boundary  which  provides  the  planetary 
scale  differential  heating.  A  number  of  solutions  were  found  for  a  balance  between 


")  Hov^ever.  Mirakami  (personal  communication)  found  that  this  horizontal  distribution  of  diver- 
gence is  m  good  agreement  with  the  1970  summer  data. 
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the  horizontal  advection  and  vertical  diffusion  of  heat  in  the  thermodynamic  energy 
equation,  and  one  of  them  allows  a  half-cycle  out  of  phase  shift  between  the  upper 
level  vorticity  and  divergence.  However,  since  the  monsoonal  differential  heating  in 
the  real  atmosphere  is  mainly  due  to  the  latent  heat  release  of  cumulus  convection 
which  is  maximum  in  the  middle-upper  troposphere,  it  is  uncertain  whether  Barcilon 
and  Cooper's  (1974)  theory  is  completely  applicable  to  the  real  atmosphere. 

Another  explanation  for  the  observed  vorticity-divergence  phase  relationship  was 
offered  by  Paegle  and  Paegle  (1976).  They  speculated  that  a  high  pressure  region 
with  strong  pressure  gradients  may  cause  the  geostrophic  vorticity  to  become  less 
than  one-half  of  the  negative  planetary  vorticity.  This  implies  that  the  balance 
equation  becomes  non-elliptic  and,  under  certain  rather  strong  assumptions,  leads  to  a 
rapid  vorticity  adjustment  with  strong  horizontal  divergence  occurring  in  the  negative 
relative  vorticity  region.  However,  the  resultant  steady-state  absolute  vorticity 
vanishes  and  this  efficiently  eliminates  the  generation  of  vorticity  by  divergence.  Such 
large  negative  values  of  relative  vorticity  are  not  observed  in  the  large-scale  monsoon 
flow.  In  fact  the  minimum  seasonal  mean  relative  vorticity  at  200  mb  during  the  1967 
summer  as  observed  by  Krishnamurti  is  ^  —3  x  10""^  sec"  ',  which  occurs  in  the 
vicinity  of  the  Tibetan  and  Mexican  highs  near  30°N.  This  magnitude  is  considerably 
smaller  than  that  of  the  Coriolis  parameter.  Thus  it  appears  unlikely  that  this  super- 
gradient  effect  is  relevant  to  the  planetary-scale  monsoon  problem. 

To  test  the  hypotheses  proposed  by  Holton  and  Colton  (1972).  Chang  and 
Penti.monti  ( 1977)  used  a  3-!evel  primitive  equation  model  with  a  4°  finite-difference 
resolution  in  both  the  zonal  and  meridional  directions.  The  zonal  mean  flow  is 
included  using  a  Newtonian  term  which  maintains  the  proper  north-south  temperature 
gradient  by  continuously  adjusting  each  temperature  toward  a  thermal  equilibrium 
determined  by  the  zonally  symmetric  climatology.  The  zonally  asymmetric  motions 
are  generated  by  a  specified  perturbation  heating  function  whose  horizontal  distribu- 
tion follows  that  of  the  2(J0  mb  seasonal  mean  divergence  observed  by  Krishna.murti 
(1971b).  It  turns  out  that  this  heating  distribution  almost  reproduces  the  observed 
mean  divergence  in  the  upper  troposphere,  which  suggests  thai  the  choice  of  the 
forcing  function  is  proper.  Two  experiments  were  carried  out  by  integrating  the  model 
in  time  until  25  days  after  it  reached  the  respective  steady  states.  Both  experiments 
include  only  weak  damping  (time  scale  >7  days)  in  the  free  atmosphere  but  the 
heating  function  was  held  steady  only  in  the  first  experiment  and  allowed  to  fluctuate 
in  time  with  an  amplitude  equal  to  one-half  of  the  time-mean  heating  in  the  second 
experiment.  The  fluctuation  period  was  set  to  10  days.  Comparing  the  results  of  the 
two  experiments,  Chang  and  Pentimonti  found  essentially  no  difference  between 
the  time-mean  planetary  scale  fields  except  a  slightly  reduced  amplitude  for  the 
fluctuating  case.  Figure  7  is  their  time-mean  wind  field  at  250  mb  ( near  the  upper  level 
of  the  model)  for  the  fluctuating  heating  case.  It  can  be  seen  that  the  simulated 
Tibetan  and  Mexican  anticyclones  are  shifted  about  ?0°  and  25°.  respectively,  from 
the  observed  positions.  These  phase  shifts  are  almost  identical  to  those  produced  by 
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their  steady-heating  case.  Thus  the  fluctuating  component  as  specified  in  the  model 
is  almost  inconsequential  in  providing  the  damping  needed  tor  the  phase  shift  prob- 
lem. On  the  other  hand,  in  both  experiments  the  amplitudes  of  the  long  waves  are  not 
excessively  large  compared  to  observations.  This  result  is  similar  to  that  of  Abbott 
(1973)  and  seems  to  suggest  that  nonlinear  processes  could  be  responsible  for  part  of 
the  planetary-scale  vorticity  sink  deduced  by  Holton  and  Colton's  linear  calculations. 
The  fluctuation  of  the  planetary-scale  flow  apparently  enhances  these  nonlinear 
conversions  resulting  in  the  further  reduction  of  its  time-mean  amplitude. 

A  noteworthy  feature  in  Chang  and  Pentim(jnti"s  ( 1977)  time-mean  fields  is  that 
the  positions  of  the  upper  tropospheric  mid-oceanic  troughs  do  not  have  significant 
phase  errors  as  those  of  the  anticyclones.  As  pointed  out  by  Krishna.murti,  Daggu- 
paty.  Fein.  Kanamitsi'  and  Lee  (1973a).  the  damping  due  to  cumulus  transports 
should  not  operate  over  the  relatively  cloud-free  oceanic  areas.  With  the  inclusion  of 
nonlinear  processes  (and  also  large-scale  vertical  advections),  the  damping  require- 
ment implied  by  the  phase  errors  of  Chang  and  Pentimonti's  time-mean  upper  level 
anticyclones  is  weaker  than  that  found  by  Holton  and  Colton  (1972).  Perhaps  a 
damping  time  longer  than  1  day,  which  seems  more  reasonable  than  Holton  and 
Colton's  value  if  interpreted  as  due  to  the  cumulus  transports,  may  be  sufficient. 
In  this  case  a  numerical  experiment  similar  to  Chang  and  Pentimonti's,  but  in- 
corporating a  moderately  strong  damping  term  with  time  scale  order  of  a  few  days  in 
the  continental  convective  area  only,  may  provide  some  insight  as  to  whether  cumulus 
transport  processes  are  responsible  for  the  missing  vorticity  sink. 

The  studies  mentioned  above  all  concern  the  summer  monsoon  only.  During 
winter  Krishnamurti  et  al.  (1973b)  observed  that  the  rotational  component  of  the 
seasonal-mean  flow  of  the  active  monsoon  region  is  almost  one  order  smaller  than 
that  of  the  summer  monsoon,  while  the  divergent  component  has  a  comparable 
magnitude.  The  reason  for  this  appears  to  be  the  proximity  of  the  winter  monsoon 
convections  to  the  equator,  which  reduces  the  elficiency  of  vorticity  generation  by 
horizontal  divergence.  However,  the  importance  of  a  dissipation  mechanism  in  the 
winter  monsoon  was  also  brought  out  in  a  few  theoretical  studies  recently.  In  a  linear 
study  of  the  stationary  forcing  in  the  tropical  atmosphere.  Webster  (1972,  1973) 
found  that  the  primary  response  in  the  equatorial  region  corresponds  to  the  atmos- 
pheric Kelvin  waves  (Holton  and  Lindzen.  1968).  A  theoretical  calculation  shows 
that  the  zonal  wavelength  of  this  Kelvin  wave  response  is  very  long  while  its  vertical 
wavelength  is  very  short  [0(1  km)].  The  former  agrees  quite  well  with  observed 
equatorial  stationary  circulations  such  as  the  winter  monsoon,  but  the  latter  is  one 
order  of  magnitude  shorter  than  the  observed  vertical  wavelengths.  This  discrepancy 
arises  because  the  equatorial  Kelvin  waves  behave  like  internal  gravity  waves  in  the 
zonal-vertical  plane  and  their  phase  speed,  r,  is  related  to  an  equivalent  depth  scale  /; 
by  the  relationship  c  =  U/h)'  '.  These  stationary  waves  are  imbedded  in  a  mean 
easterly  zonal  current  so  that  thev  have  a  Doppler-shifted  phase  speed  of  a  few  meters 
per  second,  which  implies  a  very  short  vertical  wavelength.  To  get  around  this  problem 
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in  simulating  the  observed  stationary  motions,  Webster  (1972,  1973)  adopted  a  two 
level  numerical  model  with  friction  which  is  quite  successful  in  reproducing  many  of 
the  observed  features. 

Very  recently  Chang  (1977)  re-examined  this  problem  by  incorporating  linear 
damping  terms  in  the  thermally  forced  Kelvin  wave  equations.  The  dispersive  rela- 
tionship of  the  Kelvin  waves  may  be  written  as  ck  =  5''-,  where  J.  is  the  vertical 
wavenumber  and  S  is  a  static  stability  parameter.  The  linear  damping  introduces  an 
imaginary  part  to  c  which  causes  a  to  become  complex  also.  Consequently  the 
dispersive  relationship  becomes  quadratic  and  has  two  types  of  solutions.  The  first 
type  arises  when  the  inertial  time  scale  is  faster  than  the  damping  time  scale.  This 
solution  resembles  the  regular  internal  gravity  waves  whose  vertical  wavelength 
decreases  as  the  Doppler-shifted  phase  speed  decreases.  The  second  type,  called  the 
viscous  mode  by  Chang  ( 1977),  arises  when  the  damping  time  scale  dominates  over 
the  inertial  time  scale  and  its  vertical  wavelength  increases  as  the  Doppler  phase  speed 
decreases  ( Fig.  8).  This  solution  is  also  much  more  trapped  in  the  vertical  so  that  it  has 
appreciable  amplitude  only  in  the  region  of  thermal  forcing  (troposphere).  In  the 
inviscid  limit  this  solution  reduces  to  a  trivial  solution  with  vanishing  amplitude.  For  a 
reasonable  range  of  physical  parameters,  Chang  (1977)  calculated  a  vertical  wave- 
length of  >  17  kms  for  the  stationary  response,  which  is  a  significant  improvement 
over  the  inviscid  theory.  The  horizontal  structure  of  this  solution  is  also  in  good 


Figure  8 
Real  Doppler-shifted  phase  speed,  c^.  as  a  function  of  vertical  wavelength,  L..  and  damping  tnne  for  the 

Kelvin  waves.  (From  Chang.  1977.) 
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agreement  with  observations  and  Webster's  (1972,  1973)  numerical  model.  This 
requirement  of  a  damping  mechanism  with  time  scale  on  the  order  of  5  days  is  much 
weaker  than  that  found  by  Holton  and  Colton  (1972)  for  the  summer  monsoon 
and  could  well  be  explained  by  cumulus  transports.  Detailed  observational  studies  of 
the  winter  monsoon  are  needed  to  improve  the  theoretical  understanding  of  this 
problem. 


3.  Localized  hurotropic  instuhility  of  the  planetary  scale  summer  monsoon 

As  mentioned  previously,  the  amplitude  of  stationary  long  waves  during  summer 
is  better  simulated  by  the  nonlinear  models  than  by  the  linear  models  when  excessively 
large  damping  is  excluded.  This  apparently  is  due  to  the  barotropic  energy  conversions 
from  planetary  to  both  zonal  and  synoptic  scales  which  are  allowed  in^the  nonlinear 
models  only.  Observational  evidence  of  this  type  of  energy  conversion  has  been 
presented  by  Kanamitsl'.  Krishnamurti  and  Depradine  ( 1972).  who  computed  the 
wave-zonal  ilow  and  wave-wave  energy  exchanges  for  the  tropical  belt  from  15°S  to 
15°N  using  the  1967  summer  200  mb  wind  data  analyzed  by  Krishnamurti  and 
Rogers  ( 1 970).  In  particular,  they  showed  that  the  synoptic  scale  waves  receive  kinetic 
energy  at  a  fairly  significant  rate  (order  of  10"^  m-  sec  -^)  from  both  the  zonal  flow 
and  the  planetary-scale  waves.  A  similar  result  for  the  equator  to  3(rN  belt  was  also 
obtained  by  Kri.shna.mlrti  et  al.  (1973b). 

Transient  synoptic-scale  wave  disturbances  have  been  observed  in  the  tropical 
upper  troposphere  in  all  seasons.  During  summer,  the  most  active  season,  the  distur- 
bances are  found  to  occur  primarily  in  two  types  of  regions.  One  type  is  in  the  vicinity 
of  the  mid-oceanic  troughs  (Rick.  1959:  Sadler.  1967;  Miller  and  Carlson.  1970; 
Frank.  1970)  where  the  disturbances  may  assume  the  form  of  open  waves  or  closed 
vortices.  Occasionally  these  disturbances  may  develop  in  the  vertical  or  be  coupled 
with  low  level  waves,  and  intensification  leading  to  tropical  cyclone  development 
sometimes  may  also  follow.  Another  region  of  transient  wave  activity  is  along  the 
upper  tropospheric  easterly  jet  over  the  Indian  ocean,  where  Krishna.murti  and 
Rogers  ( 1 970)  have  reported  fairly  regular  occurrences  of  synoptic-scale  disturbances. 
Apparently  these  regions  are  where  most  of  the  barotropic  energy  conversion  from 
zonal  and  planetary  scales  to  the  synoptic  scale  takes  place.  Thus  it  appears  that  the 
nonlinear  energy  transfers  of  the  summer  monsoon  have  a  persistently  localized 
nature.  Another  evidence  of  this  is  related  to  the  linear  barotropic  instability  criterion, 
which  has  been  shown  by  Merilees  (1968)  to  be  usually  satisfied  if  nonlinear  energy 
exchanges  are  very  prominent.  The  upper  tropospheric  zonal  mean  flow  of  the 
summer  monsoon  does  not  satisfy  this  criterion,  but  regions  of  large  negative 
gradients  of  relative  vorticity  (Fig.  9)  where  this  instability  criterion  is  satisfied 
locally  can  be  found  in  the  seasonal  mean  200  mb  analysis  (Krishnamurti.  1971b). 
These  negative  vorticity  gradients  occur  mainly  in  the  vicinity  of  the  easterly  jet  and 
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Figure  9 
Summer  mean  (June-August  1967)  relative  vorticity  at  200  mb.  contour  inter.al  is  10"''  sec"'.  (From 

Krishnamurti,  1971b.) 

the  mid-oceanic  troughs,  and  on  some  individual  days  their  magnitude  may  even 
greatly  exceed  the  seasonal  mean  values.  Therefore  it  is  plausible  that  at  least  some  of 
the  transient  synoptic  wave  activities  in  these  regions  is  generated  by  the  localized 
(and  sometimes  short  term)  barotropic  instability,  which  is  clearly  a  result  of  the 
zonal  asymmetry  of  the  planetary-scale  summer  monsoon. 

The  kinetic  energy  conversions  from  the  planetary  scale  and  the  zonal  flow  to  the 
synoptic  scale  in  the  upper  troposphere  have  been  simulated  numerically  by  Abbott 
(1973)  and  Krishnamurti  ei  al.  (1973a).  The  conversions  produced  by  Abbott's 
(1973)  model  are  roughly  of  the  same  order  of  magnitude  as  thos?  computed  by 
Kanamitsu  el  al.  ( 1 972).  However,  the  conversions  obtained  by  Krishnamurti  et  al. 
(1973a)  are  about  one  order  smaller.  This  appears  to  be  related  to  their  treatment  of 
the  Tibetan  high  as  a  fixed  barrier  which  prevents  vorticity  transport  across  the 
boundary  of  the  barrier.  As  a  consequence,  barotropic  instability  cannot  be  realized 
in  the  immediate  vicinity  of  the  Tibetan  high  and  the  total  amount  of  energy  extracted 
from  the  larger  scales  by  the  synoptic  scale  within  a  fixed  period  of  time  must  be 
greatly  reduced.  Nevertheless,  Krishna.murti  et  al.  (1973a)  were  able  to  simulate  the 
gross  features  of  the  planetary  scale  waves  in  the  upper  troposphere,  such  as  the 
mid-oceanic  troughs  and  the  Mexican  high.  Apparently  once  the  strength,  shape  and 
position  of  the  Tibetan  high  are  fixed  in  the  experiment,  the  interaction  of  the  zonal 
mean  flow  and  this  blocking  high  is  sufficient  to  produce  mechanically  many  other 
features  of  the  circulation. 

In  the  numerical  model  by  Colton  ( 1973),  attention  is  focused  on  the  development 
of  the  synoptic-scale  disturbances  in  the  easterly  jet  and  in  the  mid-Pacific  trough 
regions  where  the  barotropic  instability  criterion  is  locally  satisfied.  Near  the  easterly 
jet  the  model  produces  disturbances  that  propagate  westward  at  a  phase  speed  slightly 
less  than  the  local  speed  of  the  jet  and  continue  to  grow  until  they  are  well  past  the 
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longitude  of  the  jet  maximum.  After  reaching  the  maximum  amplitude  they  begin  to 
decay  fairly  rapidly.  Many  features  of  these  disturbances  resemble  those  analyzed  by 
Krishnamurti  and  Rogers  (1970)  in  the  Indian  Ocean  and  African  region.  In  the 
mid-Pacific  trough  the  model  also  produces  waves  and  vortices  whose  movement  and 
behavior  are  quite  similar  to  those  observed  in  daily  synoptic  charts.  In  both  regions 
the  growing  disturbances  exhibit  a  latitudinal  tilt  in  opposition  to  the  lateral  shear  of 
the  quasi-steady  basic  flow,  which  is  characteristic  of  the  barotropically  unstable 
waves. 


I 


Figure  10 

Zonal  average  kinetic  energy  at  20N  as  a  function  of  time  and  wavenumber  as  computed  by  Colton(  1973). 

Contour  interval  is  5  m'  sec"".  Hatched  areas  are  regions  of  energy  maximum. 
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In  studying  the  barotropic  energy  exchanges.  Colton  (1973)  also  examined  the 
time  evolution  of  the  computed  kinetic  energy  spectrum.  An  example  of  this  is  given 
in  Fig.  10  which  shows  the  zonally  averaged  energy  at  20°N  as  a  function  of  time  and 
wavenumber.  A  noteworthy  feature  is  that  when  the  energy  of  wavenumber  1  is  at  a 
maximum,  the  energy  at  the  short  waves  is  usually  at  a  minimum,  and  vice  versa.  This 
suggests  that  energy  exchanges  take  place  continuously  between  the  forced  long  waves 
and  the  transient  disturbances.  The  disturbances  remove  energy  from  the  planetary 
scale  during  the  growing  stage  and  feed  energy  back  into  the  planetary  scale  during  the 
decaying  stage.  As  a  consequence  the  amplitude  of  the  time  mean  planetary-scale 
waves  is  smaller  than  that  produced  by  a  linear  model.  Since  no  dissipation  is  provided 
for  the  transient  waves  all  energy  extracted  from  the  forced  motion  must  either 
remain  in  the  transient  waves  or  be  transferred  back  to  the  forced  motion.  Colton 


bL 


I 


1102  Chih-Pei  Chang  (Pageoph, 

(1973)  speculated  that  the  actual  damping  of  the  planetary  scale  by  the  local  baro- 
tropic  instability  may  be  more  effective  if  some  dissipation  mechanism  for  the 
transient  synoptic  waves  is  included. 

Although  Colton's  (1973)  model  lacks  a  sufficient  damping  mechanism  and 
therefore  produces  a  basic  flow  that  is  somewhat  too  intense  and  has  a  phase  shift 
problem,  it  nevertheless  illustrates  the  interesting  phenomenon  of  waves  developing 
and  decaying  in  a  spatially  varying  basic  flow  which  contains  locally  unstable  and 
stable  regimes.  To  gain  further  physical  insight  into  this  phenomenon,  especially  with 
respect  to  the  easterly  jet  of  the  summer  monsoon,  Tupaz,  Williams  and  Chang 
(1977)  formulated  a  1-level  numerical  model  using  the  barotropic  vorticity  equation 
to  study  this  problem  from  a  linear  instability  point  of  view.  In  their  model  the  basic 
zonal  wind  is  approximated  by  a  steady  easterly  Bickley  (hyperbolic-secant-squared) 
jet  which  increases  slowly  from  the  inflow  (eastern)  boundary  to  a  longitude  in  the 
middle  of  a  rectangular  domain,  then  decreases  slowly  downstream  towards  the 
outflow  (western)  boundary.  The  basic  meridional  wind  is  derived  in  such  a  way  that 
the  total  basic  flow  is  non-divergent.  Periodic  forcing  is  introduced  on  the  inflow 
boundary  and  a  radiation  condition,  whereby  a  wave  is  allowed  to  propagate  its 
energy  out  of  the  domain,  is  applied  at  the  outflow  boundary.  The  model  is  integrated 
from  zero  interior  perturbation  streamfunction  until  the  perturbation  fields  vary 
periodically  in  time  at  each  grid  point.  At  this  quasi-steady  state  the  waves  move 
through  the  region  and  grow  or  decay  in  space  only,  in  reaction  to  the  local  stability 
properties  of  the  basic  flow. 

To  compare  the  numerical  results  with  the  "local  growth  rate"  concept  of  the 
parallel  flow  theory,  Tupaz  et  al.  (1977)  developed  a  simple  analytical  model  to 
calculate  the  effects  of  the  local  instabilities.  This  is  done  by  considering  an  equation 
that  allows  for  propagation  and  growth  or  decay,  i//,'  +  c{x)il/'^  =  /j(  v)i//'.  where  /  is 
time.  X  the  zonal  coordinate,  and  i//'  the  perturbation  streamfunction.  The  local 
phase  speed  c  and  local  growth  rate  //  can  be  obtained  from  the  parallel  flow  theory. 
When  a  periodic  solution  with  a  frequency  (o  and  an  amplitude  F(x)  is  introduced,  the 
downstream  variation  of  f  with  respect  to  a  reference  point  x^  may  be  written  as 

F{x)  =  F(.Vo)exp[-/        w  c/.v/c(.v)]  exp  [       n{x)  dx/c{x)1. 

J  Xll  J-Vo 

In  this  equation  the  local  wavelength  is  co/c  and  the  local  spatial  growth  rate  is  n/c. 
Figure  1 1  compares  the  steady  state  envelope  for  i//'  of  the  numerical  solution  and  the 
corresponding  F(x)  for  a  case  in  which  the  central  speed  of  the  easterly  jet  is  increased 
from  15  m  sec" '  at  V  =  —  10  (X)0  km  to  30  m  sec"  '  at  v  =  0  and  then  reduced  to 
15  m  sec"  '  at  v  =  10  000  km.  The  curve  for  F(x)  can  be  shifted  up  or  down  in  the 
diagram  by  changing  Xq,  but  its  shape  does  not  change.  Two  F{x)  curves  are  included 
in  Fig.  1 1.  one  for  .Xo  =  280  km  and  the  other  for  Xq  =  4200  km,  both  based  on  the 
local  wavelength  distribution  measured  from  the  numerical  solution.  It  is  evident 
from  this  figure  that  both  the  numerical  solution  and  the  local  growth  rate  of  the 
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Figure  1 1 
Comparison  of  wave  packet  envelopes  in  a  locally  unstable  mean  flow:  ' li/    ■.  computed  from  a  numerical 
model ;  f ,(  v„  =  2X0  km  i  and  F,,l  v„  =  4200  km),  computed  from  an  analytical  model  based  on  the  parallel 

flow  theory.  (From  Tupaz  et  al.  19''?.) 

parallel  flow  theory  put  the  ma.ximum  amplitude  of  the  perturbations  at  .v  =  —9800 
km.  well  past  the  jet  maximum.  This  is  in  good  agreement  with  the  results  obtained  by 
CoLTON  (1973)  and  Chang  and  Pentimonti  (1977).  and  is  qualitatively  consistent 
with  KRiSHNAMURTiand  Rogers"  ( 1970)  observation  that  the  most  active  disturbance 
area  extends  downstream  from  the  jet  maximum  near  the  southern  tip  of  India  to 
central  Africa  (see  Fig.  10.2,  top  panel,  o^  Krishna.vu'RTI,  1971b).  This  longitude  of 
maximum  amplitude  corresponds  to  the  neutral  point  of  the  basic  flow,  and  is 
immediately  understandable  from  the  local  growth  rate  concept:  a  wave  coming  in 
from  the  inflow  boundary  will  grow  as  it  moves  downstream  in  the  locally  unstable 
region,  until  it  reaches  the  neutral  point  where  it  begins  to  decay.  However,  the 
numerical  solution  has  a  larger  maximum  amplitude  than  that  predicted  by  the  local 
growth  rate  of  the  parallel  flow  theory.  The  waves  actually  grow  faster  in  the  unstable 
region  and  decay  faster  in  the  stable  region.  Ti'Paz  et  al.  (1977)  evaluated  the  energy 
equation  for  the  numerical  model  and  found  that  the  most  important  term  for 
perturbation  energy  production  is  the  Reynold  stress  term  proportional  to  the 
meridional  shear  of  the  basic  zonal  wind.  This  term  is  also  the  only  source  term  for  a 
parallel  flow  model,  and  it  depends  on  the  phase  structure  of  the  perturbation  held. 
Figure  12  compares  the  local  phase  tilt  of  the  wave  from  the  numerical  model  with 
that  from  the  parallel  flow  theory  at  two  longitudes:  v  =  0  km  (jet  maximum)  and 
.v  =  -2800  km  (slightly  downstream  from  jet  maximum).  At  both  places  the  phase 
tilt  from  the  numerical  model  is  significantly  larger  than  that  of  the  parallel  flow  theory. 
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Thus  apparently  the  downstream  variation  of  the  basic  flow  augments  the  phase  tilt 
of  the  waves  which  gives  a  larger  growth  rate  in  the  unstable  region.  A  similar  situation 
is  found  in  the  stable  region  where  the  dynamic  damping  due  to  the  continuous 
spectrum  effect  is  enhanced  by  a  stronger  lilt  in  the  same  direction  of  the  basic  shear  of 
the  easterly  jet.  These  results  indicate  that  the  barotropic  energy  exchanges  associated 
with  a  zonally  varying  easterly  jet  which  is  locally  unstable  can  actually  be  more 
prominent  than  those  predicted  by  the  local  growth  rate  from  the  parallel  flow  theory. 
Hence  the  presence  of  the  planetary-scale  asymmetries  during  the  summer  monsoon 
strengthens  the  nonlinear  energy  conversions  resulting  in  significant  activity  of 
synoptic-scale  disturbances. 
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Figure  12 

Comparison  of  the  phase  tilts  of  waves  in  a  locally  unstable  mean  flow  computed  by  a  numerical  model 

(8*)  and  the  parallel  flow  theory  (('„)  at  two  longitudes:  (a)  x  =  0  where  u{x.  0)  =  -30  m  sec    '  and 

(b)  .X  =  28(X)km  where  u(x.  0)  =   -27.3  m  sec''.  (From  Tupaz  f/ u/.  1977.) 


Although  TuPAZ  et  al.  (1977)  employed  linear  equations  only,  their  results  may 
also  shed  some  light  on  the  eff"ects  of  the  transient  waves  on  the  basic  flow.  The  waves 
remove  kinetic  energy  from  the  basic  flow  (both  zonal  and  planetary  scale)  and  most 
of  this  energy  is  removed  on  the  downwind  side  of  the  jet  ma.ximum.  Therefore  it  may 
be  expected  that  in  a  nonlinear  extension  of  their  model  the  planetary  scale  will 
experience  maximum  damping  west  of  the  easterly  jet  maximum,  which  undoubtedly 
will  affect  the  amplitude  and  phase  of  the  forced  planetary-scale  waves. 
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4.   Oscillations  of  the  planetary-scale  monsoons 

The  summer  monsoon  circulation  is  known  to  exhibit  many  fluctuations  in  a 
variety  of  time  scales  within  each  season.  The  most  notable  fluctuation  is  the  alter- 
nating active  and  break  periods  of  the  monsoon.  During  the  break  periods  most  of 
the  circulation  characteristics  may  diminish  or  even  reverse  over  the  broad  monsoon 
region.  Less  dramatic  variations  occur  more  frequently  than  the  active-break  cycles, 
although  sometimes  they  are  limited  to  local  regions.  Figure  13  shows  the  power 
spectrum  of  the  1967  summer  200  mb  streamfunction  over  a  region  of  the  Tibetan 
high,  analyzed  by  Krishnamurti  et  al.  ( 1973a).  A  pronounced  oscillation  periodicity 
is  noted  around  10-13  days  and  a  secondary  spectral  peak  occurs  near  3  days.  A 
similar  spectral  distribution  for  the  tropospheric  wind  field  of  the  1962  summer  was 
found  by  M.  Muraka.mi  (1976),  where  the  main  peaks  occur  at  10-15  days  and  4-5 
days.  Both  of  these  investigators  felt  that  the  shorter  period  oscillation  may  be  a 
manifestation  of  the  synoptic-scale  westward  propagating  monsoon  depressions.  The 
longer  period  oscillation,  on  the  other  hand,  appears  to  be  related  to  the  break  monsoon 
and  other  broad-scale  variations  of  the  monsoon  circulation.  The  10-15  day  periodi- 
city of  this  oscillation  has  also  been  noted  during  many  other  years.  For  example,  it  is 
reflected  in  the  daily  rainfall  data  of  1963.  1971  (Murakami.  1972),  and  1972,  1973 
(Krishnamlrti  and  Bhalme.  1976). 

After  detecting  a  spectral  peak  near  the  13-day  periodicity  in  many  parameters  of 
the  planetary-scale  monsoon  system  of  the  1967  summer,  Krishnamurti  and  Bhalme 
(1976)  calculated  the  phase  relationship  between  the  oscillatioji  of  these  parameters. 
Their  result  suggests  that  in  a  13-day  cycle  the  sequential  strengthening  of  the  mon- 
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soon  system  including  the  monsoon  trough,  the  cloud  brightness,  the  surface  Mas- 
carene  high,  the  upper  Tibetan  high  and  easterly  jet,  and  the  central  Indian  rainfall 
occur  after  the  dry  and  moist  destabilization  of  the  lower  troposphere.  This  will  then 
be  followed  by  a  weakening  of  the  system  in  the  same  sequence  after  dry  and  moist 
stabilization.  On  the  other  hand,  M.  Murakami  ( 1976)  used  a  composite  technique 
to  study  the  1962  summer  monsoon,  noting  that  it  is  difficult  to  verify  the  statistical 
significance  of  the  cross-spectral  quantities  of  this  period  range  derived  from  a 
3-month  time  series  of  daily  data.  He  found  that  the  intensification  and  weakening  of 
convection  and  other  low-  and  upper-level  winds  are  closely  in  phase  for  the  1 5-day 
cycle  during  this  season.  In  any  case,  there  is  strong  evidence  that  many  components  of 
the  planetary-scale  monsoon  are  involved  in  this  oscillation. 

The  observed  10-15  day  oscillation  raises  an  interesting  question  concerning  its 
physical  mechanism.  Krishnamurti  and  Bhalme  ( 1976)  proposed  that  it  may  be  a 
result  of  a  cloud-radiative  feedback  mechanism  in  which  the  enhanced  convection 
causes  increased  cloud  cover  which  in  turn  reduces  the  incoming  short-wave  radiation 
and  stabilizes  the  lower  troposphere.  This  is  then  followed  by  decreased  convection 
and  cloudiness,  and  therefore  increased  short-wave  radiative  heating  at  the  surface 
which  renews  the  active  convection.  Another  possible  oscillation  mechanism  is 
revealed  by  a  simple,  domain-averaged,  atmosphere-ocean  interacting  model  by 
Webster  and  Lau  ( 1977).  In  their  model  the  large  heat  capacity  of  the  ocean  causes 
it  to  respond  slowly  lo  the  solar  heating  while  the  atmosphere  temperature  reaches  a 
ma.ximum  due  to  heating  at  the  land  surface.  After  a  period  of  time  the  ocean  tem- 
perature peaks  and  this  feeds  back  to  the  atmosphere  through  convection  resultmg  m 
another  temperature  ma.ximum.  Although  this  process  alone  is  unlikely  to  produce 
subsequent  maxima  without  the  participation  of  other  mechanisms,  it  may  still  be 
relevant  to  some  of  the  active-break  cycles  because  Rama.mlrth'^'s  (1969)  survey 
(quoted  in  Krishnamurti  and  Bhal.me,  1976)  indicates  that  there  are  only  one  or  two 
major  breaks  during  most  oi'  the  summer  monsoon  seasons.  Other  possible  mechan- 
isms for  the  planetary-scale  oscillation  may  include  interactions  between  hemispheres 
and  between  stationary  and  propagating  long  waves  as  mentioned  by  Krishnamurti 
and  Bhalme  (1976),  and  nonlinear  vacillations  between  different  scales.  The  latter 
mechanism  produces  a  5-day  oscillation  of  the  planetary  scale  in  Colton's  (1973) 
barotropic  model,  which  is  purely  dynamical  and  does  not  provide  dissipation 
processes  for  the  synoptic  scale.  If  thermodynamics  and  synoptic  damping  are 
included,  the  interaction  would  become  more  complex  and  the  oscillation  period 
could  be  modified. 

The  planetary-scale  circulation  during  the  winter  monsoon  also  exhibits  significant 
variations,  largely  due  to  the  influence  of  the  cold  air  outbreaks  from  the  Asian 
continent.  These  outbreaks  often  interfere  with  the  main  winter  con\ective  region  in 
the  maritime  continent  and  cause  widespread  weather  changes  in  the  South  China 
Sea  and  Vlalaysia-lndonesia  region,  including  the  development  of  strong  synoptic- 
scale  disturbance  activity.  Since  this  convective  region  provides  the  major  energy 
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source  for  the  thermally  direct,  planetary-scale  east-west  overturning  and  Hadley 
cells  (Krishnamurti  et  al.  1973b),  its  variations  directly  influence  the  circulations  in 
the  entire  tropics  and  southern  mid-latitudes  during  winter.  These  planetary-scale 
zonal  and  lateral  teleconnections  are  reflected  in  the  analyses  by  Murakami  and 
Unninayar  (1977)  and  Krishnamurti  et  al.  (1973b)  despite  their  limited  data  base. 
The  initiation  and  progress  of  theoretical  studies  in  this  area  would  depend  heavily  on 
the  expansion  of  our  present  observational  knowledge  of  this  problem. 


5.   Concluding  remarks 

The  three  problems  reviewed  here  are  inter-related  to  each  other.  For  example, 
the  localized  barotropic  instability  is  a  reflection  of  the  strong  nonlinear  energy 
transfer  which  appears  to  be  quite  relevant  to  the  damping  of  the  planetary-scale 
flow  during  the  summer  monsoon.  Also  Chanc;  and  PentimontTs  ( 1977)  experiment 
indicates  that  these  energy  conversions  are  more  pronounced  when  the  planetary 
scale  contains  an  appreciable  oscillation.  The  enhanced  upper  tropospheric  in- 
stability in  summer  due  to  both  time  and  spatial  variations  of  the  planetary  scale  also 
has  important  ramifications  for  the  synoptic-scale  disturbances.  At  least  in  the  mid- 
oceanic  trough  regions,  the  upper-level  disturbances  may  develop  in  the  vertical  and 
could  be  related  to  some  o'i  the  weather-producing  easterly  waves  and  tropical 
cyclones.  A  different  type  of  problem  arises  for  the  winter  monsoon  during  which  the 
upper  tropospheric  synoptic-scale  waves  and  zonal  flow  feed  energy  to  the  planetary 
scale  (Krishnamurti  ct  al.  1973b).  This  on  the  one  hand  removes  the  nonlinear 
processes  as  a  possible  damping  mechanism,  and  on  the  other  hand  requires  an 
explanation  for  the  energy  source  of  the  synoptic  scale.  Clearly  much  more  effort  is 
needed  to  answer  these  questions,  whose  importance  lies  in  the  fact  that  a  thorough 
understanding  of  the  planetary  scale  is  crucial  for  the  proper  simulation  of  the  entire 
multiple-scale  monsoon  circulation. 
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